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MAPPING A TUNGSTEN -REACTOR ROCKET ENGINE AS A 
GUIDE TO OPERATION AND CONTROL 
by Richard H. Cavicchi 
Lewis Research Center 

SUMMARY 

A nonlinear analog simulation of a preliminary tungsten water-moderated reactor 
system indicated stability without external feedback controls. The simulation was used 
to obtain steady- state operating maps over a mass-flow range from 20 to 120 percent. 

A feasible operating corridor, outlined by locating system constraints on the maps, is 
bounded by the bleed-pump turbine speed limit at about 118 percent power and by bleed- 
pump stall at 33 percent power. Within this corridor, system operation is confined to a 
500° R range in chamber temperature by heat-exchanger icing and the nozzle-wall tem- 
perature limit. 

The manner of operating the bleed system determined the locations of the constraints 
on the maps except for the nozzle-wall and fuel -temperature limits. As chamber temper- 
ature was increased, temperature limits in the nozzle -wall, fuel, water-pump turbine - 
inlet, and pressure-tube were reached successively. Design specific impulse was poten- 
tially obtainable above 37 percent power, but the margin from the nozzle-wall tempera- 
ture limit was small. 

Reactor power could be controlled by the topping-turbine bypass valve and poison re- 
activity. Also, water temperature was a suitable control parameter. Reactor power 
was approximately proportional to hydrogen flow rate. A 10-cent insertion of poison re- 
activity reduced reactor power by about 7 percent; whereas less than a 3 percent change 
in power was accompanied by a 1° R change in water temperature. 

The water-pump turbine valves can be replaced by fixed orifices sized for 100- 
percent design-point operation. Further simplification in control is possible above 
47 percent power by keeping the bleed-pump turbine valves at their 100 percent settings. 
Below 47 percent power, the bleed-pump turbine valves must be controlled to delay the 
onset of pump stall. 

The bleed-gas flow rate required no control in the range of this study. This system 
revealed a self -regulating tendency in which percent bleed rose as power and total hydro- 
gen flow were decreased. In a modified system that bypassed some of the hydrogen 
around the heat exchanger, the operating envelope enlarged considerably and feasible 
steady-state operation was possible down to 13 percent power. Transient operation below 
this level may permit startup. 
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INTRODUCTION 


Extensive theoretical studies such as reference 1 have indicated the advantages of 
nuclear-powered vehicles in interplanetary missions. References 1 and 2 conclude that* 
the weight required in orbit about Earth for manned Mars missions can be at least halved 
if nuclear propulsion is used instead of chemical propulsion. 

Recent successful experimental reactor tests in the Nuclear Engine for Rocket Vehi- 
cle Application (NERVA) program (ref. 2) have increased the technical feasibility of nu- 
clear propulsion in space. The NERVA reactor is a homogeneous, graphite -moderated 
reactor. Of the many problems surmounted in the NERVA program, some of the more 
troublesome were maintaining structural integrity of the core and combating the corro- 
sive effects of chemical interaction of the hydrogen propellant with the graphite moder- 
ator (ref. 2). 

In further support of the nuclear rocket effort, NASA has undertaken a thorough in- 
vestigation of the heterogeneous class of reactor. The particular design receiving much 
attention is a thermal reactor having fuel elements of uranium mixed with tungsten en- 
riched with tungsten- 184 and using circulating water as the moderator. A feature of this 
concept is the negative temperature reactivities of the fuel and water moderator. This 
self -regulating characteristic not only acts as a built-in safety factor, but should aid in 
reactor-power control in response to any prescribed demand. Reference 3 points out 
other features, among which are low weight; high -temperature thermal-resistant charac- 
teristics and low neutron-absorption of the tungsten- 184; and relief of the moderator 
structural problems by using a fluid moderator. The corrosion problem of interaction 
between propellant and moderator becomes nonexistent. Meanwhile, the integrity of the 
fuel-element design under extremes of temperature, radiation, and aerodynamic loading 
was investigated experimentally (ref. 4). 

Aside from such reactor technical problems, the questions of stability, control, and 
operating feasibility of the integrated nuclear -engine system remained to be answered. 
For this purpose, NASA has made a complete system design with the tungsten water- 
moderated reactor as a power source in a combination bleed and topping cycle. In addi- 
tion to the reactor, the system includes a water -hydrogen heat exchanger; a bleed sys- 
tem; bleed-driven water, hydrogen, and poison turbopumps; a topping turbopump; a re- 
generatively cooled thrust nozzle; a radiation shield; circulating poison; and associated 
valves and piping. 

A preliminary version of the system design was used in the evolution of the program 
to improve the design and performance of the system. The system studied in the present 
investigation differs from a later version only by the absence of a cold-hydrogen gas pas- 
sage around the heat exchanger. 

As part of this program, the preliminary reference system was simulated on analog 
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computers so that the operation and control characteristics and interactions among the 
system parameters could be studied. This report presents the techniques used in the 
simulation and evaluates the accuracy of the results. A comparison of the computer re- 
sults with detailed digital and hand calculations made individually of the heat exchanger, 
core, and turbopumps was employed in the validation. 

The simulation was used in the present investigation to map system performance 
over a wide range of steady-state operation. By interpreting the steady-state maps ob- 
tained, the present report draws conclusions concerning the following aspects of system 
operation and control: 

(1) Location and range of the feasible operating region 

(2) Conditions restricting the operating range 

(3) Search for suitable control parameters 

(4) Merits and sensitivity of possible control parameters 

(5) Simplifications in control by fixing system variables where possible 

(6) Interactions among system parameters 

(7) Suggested design modifications for improvement in system performance 

The present analysis is one -dimensional, using a nonlinear analog computer simula- 
tion scaled primarily for use in the reactor power and flow range above 20 percent. 
Steady-state mapping was conducted by using hydrogen mass-flow rate and chamber tem- 
perature as independent variables. These parameters were varied over a range from 
20 to 120 percent and from 67 to 135 percent of design, respectively. The wide range of 
variables gave assurance of encountering all constraints and even successive limits. The 
present study does not include startup, shutdown, malfunction, or transient operation. 


DESCRIPTION OF SYSTEM 

A brief description of the preliminary system follows. A schematic diagram of the 
flow system is given in figure 1 to aid in visualizing the hydrogen and water flow paths 
and the functions of the various components. 

The heart of the system is the reactor core, which is a cylindrical aluminum tank 
containing a number of axially mounted aluminum tubes surrounded by water. Inside the 
tubes are fuel elements consisting of several concentric rings of fueled metallic cylinders 
mounted within a tungsten support tube. Hydrogen gas is heated as it flows through the 
annular passages between the rings. An aluminum pressure tube contains the concentric 
fuel-element assembly, and an annular region of stagnant hydrogen is provided between 
the tungsten support tube and the inner surface of the pressure tube. The stagnant hydro- 
gen functions as an insulator between the high-temperature hydrogen and the moderator 
water. Another aluminum tube is mounted outside the pressure tube to act as a flow 
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Figure 1. - Schematic flow diagram of tungsten water-moderated reactor preliminary reference system. 
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divider for the water. The moderator water is thus divided into a high-velocity region 
between the pressure tube and the flow divider, to enhance heat transfer, and a low- 
velocity region outside the flow divider. Several of the fuel assemblies are allocated 
£or heating hydrogen for use in driving the bleed turbines. These fuel assemblies and 
the turbines along with the associated valves are referred to as the bleed system. 


Hydrogen Flow Path 

Liquid hydrogen from the tank is fed to the first- and second- stage hydrogen pumps 
in series (see fig. 1). They are referred to herein as bleed pump and topping pump, re- 
spectively. High-pressure hydrogen passes through nozzle cooling tubes for regenerative 
cooling and on to the heat exchanger. This is a shell-and-tube parallel-flow exchanger 
mounted just inside the reactor pressure vessel. Cold hydrogen flows through the tubes 
to cool the water flowing axially through the shell. The hydrogen is collected in a plenum 
at the heat-exchanger exit and ducted to the topping turbine, which drives the second- 
stage hydrogen pump. From the topping turbine the hydrogen reenters the pressure ves- 
sel, passes through a radiation shield, and into a plenum for distribution to the fuel ele- 
ments. While flowing through the annular spaces between the fuel rings, the hydrogen 
gas is heated to an elevated temperature for discharge through the nozzle to produce 
thrust. 


Water Flow Loop 

The primary role of the circulating water is to function as the moderator in the re- 
actor core, where it thermalizes the neutrons produced by fission. In addition, it cools 
the pressure tube, flow divider, and reflector, and delivers the heat absorbed in these 
components to the cold hydrogen in the heat exchanger. 

The water loop is indicated by the dotted lines in figure 1. Water from the bleed- 
driven water pump enters an annular manifold inside the reactor pressure vessel where 
it is directed to the following three parallel paths: (1) high-velocity regions between the 
pressure tubes and the flow dividers, (2) low-velocity regions outside the flow-divider 
tubes, and (3) axial passages in the side reflector. The three water streams are col- 
lected and mixed in a plenum at the core exit. The direction of water flow is then re- 
versed as the water enters the heat exchanger and flows axially through the shell side. 
Finally, it is returned to the water pump. The water system pressure is designed to 
match the core -gas exit pressure to minimize the loads on the outlet tube sheet. 
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Bleed Turbines 


Bleed gas, heated to a suitable temperature level in the core, furnishes the energy 
to drive the two bleed turbines, which are in series. As shown in figure 1, the water- i 
pump turbine is positioned upstream of the hydrogen bleed-pump turbine, which exhausts 
to the bleed nozzle. The order of the two bleed turbines in the series arrangement was 
dictated by a low limit of practical turbine blade height; that is, if the bleed-pump tur- 
bine were situated first in the series arrangement, its high rotational-speed requirement 
in combination with the high bleed-gas pressure would demand blade heights too short to 
be practical. The lower speed requirement of the water-pump turbine, however, per- 
mits it to be placed first without needing abnormally short blades. 

During development of the turbine simulation in the present model, a preliminary 
study was made comparing parallel and series turbine arrangements. Strong interac- 
tions between the turbines were observed in response to valve changes in the series con- 
figuration. A vast improvement in turbine control became available by using a parallel 
arrangement, but this arrangement required the bleed-nozzle weight to be doubled and 
the bleed-flow requirements to be increased over 38 percent. 


Controls 

For control, the following six manipulated variables were provided in the system de- 
sign to permit maximum flexibility in operation: 

(1) Poison reactivity 

(2) Water-pump-turbine throttle valve (WTTV) 

(3) Water-pump-turbine bypass valve (WTBV) 

(4) Bleed-pump-turbine throttle valve (BTTV) 

(5) Bleed-pump-turbine bypass valve (BTBV) 

(6) Topping-turbine bypass valve (TTBV) 


SYSTEM MODEL 

The present study used a nonlinear model of the preliminary reference system in 
which the following components were simulated on analog computers: 

(1) Heat exchanger 

(2) Core 

(3) Reflector 

(4) Bleed system 
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(5) Nozzle 

* (6) Turbopumps 
(7) Valves 

* (8) Reactor kinetics 

Three E.A.I. 231-R analog computers were used in the simulation using the follow- 
ing equipment: 287 amplifiers (including 53 integrators), 85 multipliers, 27 fixed and 
variable diode function generators, and 330 potentiometers. 

The analysis, including the symbols, assumptions, basic equations, and application, 
is presented in appendixes A to D, respectively. 

TABLE I. - NUMBER OF SECTIONS USED 
IN MODEL COMPONENTS 


Component 

Equation 

Energy 

Momentum 

Energy 

Number of sections 

Axial 

Radial 

Core fuel assembly 

3 

— 

1 

Core gas 

3 

1 

1 

Core water 

1 

— 

1 

Reflector material 

1 

--- 

1 

Reflector water 

1 

--- 

1 

Bleed fuel assembly 

1 

— - 

1 

Bleed gas 

1 

1 

1 

Heat -exchanger material 

3 

— 

2 

Heat -exchanger gas 

3 

1 

1 

Heat -exchanger water 

3 

— 

1 

Nozzle coolant 

1 

1 

1 


A summary of the number of sections used in the component simulations, where ap- 
plicable, is presented in table I. As used herein, a section is defined as an imaginary 
subdivision in the representation of a component for use in applying finite-difference ap- 
proximations to differential equations. Sections may be formed in the axial or radial di- 
rections, or both. 

The simulation was scaled for the power and flow range above 20 percent. Startup 
and shutdown operations, which are not included in this study, would require rescaling 
throughout the system. 
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Heat Exchanger 


The heat exchanger was represented in the simulation by a single tube, containing 
:old hydrogen gas flow, on the outside of which hot moderator water flowed. Both fluids' 
flow in the same direction in the heat exchanger. 

Gas. - Preliminary transient analog runs indicated that at least three axial sections 
were required for the heat -exchanger gas energy equation and one for the momentum 
equation. Dynamic traces of one- and two-section energy-equation models differed 
noticeably from those obtained by using three sections. Use of two or three sections to 
simulate the momentum equation, however, yielded no improvement over the use of one 
section. The length of each of the three sections for the energy equation is one -third of 
the heat-exchanger length. 

Tubes . - Three axial sections were also used in the heat-exchanger material analy- 
sis to correspond with the three gas sections. Because of the high temperature gradient 
across the tube material (from hot water to cold hydrogen), it was found that one radial 
section was inadequate. Accordingly, two radial sections were used and the center tem- 
perature in each was calculated from the energy equation. With the assumption of a lin- 
ear variation between the two radial- section center temperatures, the tube -surface tem- 
peratures were evaluated by extrapolation. 

Water . - Three axial sections were used to simulate the heat -exchanger water 
energy equation. They correspond to the three axial sections used in the gas and tube 
material. 


Core 

Fuel assembly . - The concentric fuel rings were simulated by a single tube with a 
hydraulic diameter equivalent to that of the concentric assembly. On the basis of pre- 
liminary transient analog runs, the use of three axial sections was indicated for the rea- 
sons cited in the heat-exchanger discussion. The lengths of the first two sections are 
each one -fourth of the total length, whereas the third section is one-half the total. This 
selection of axial -section length was based on the axial reactor power distribution. One 
radial section was used. 

Gas . - Three axial sections, corresponding to those used in the fuel assembly simu- 
lation, were used for the core -gas energy representation. Preliminary transient analog 
runs revealed that one axial section was adequate to simulate the momentum equation. 

Water . - The water moderator flows through the core in units of two parallel pas- 
sages separated by tubular flow dividers of aluminum. The array of pressure-tube - 
flow-divider units, with associated parallel water streams, that enclose the individual 
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fuel elements was represented by a single unit. Transient exploratory studies on the 
analog computer indicated a preference in using a one -axial- section, two-passage model 
instead of a three -axial- section, one-passage model. The flow split is assumed to be 
constant at the 100-percent design value. 

Pressure tube and flow divider . - The pressure tube and flow divider were each 
simulated by one axial and one radial section. Therefore, the temperature of each was 
assumed constant throughout the material. 


Reflector 

The side reflector was represented by a single equivalent beryllium cylinder. A 
single passage through the reflector material represented the many passages in the 
actual design. 

Material . - One axial section was used in the reflector -material simulation. 

Water . - A one -axial- section model was also used for the reflector water. The 
water mass-flow rate through the reflector was assumed to be proportional to the total 
circulating -water flow rate. 


Bleed System 

The bleed gas makes two passes through the core. In the first pass, both the pres- 
sure drop and heat absorbed by the bleed gas are relatively small. Therefore, only the 
second pass was simulated. The shield was not simulated in detail. The heat given up 
by the bleed gas to the shield, water, and cold hydrogen at the inlet -hydrogen plenum 
was assumed proportional to the bleed-gas exit temperature. 

Fuel assembly . - A single tube was used to represent the fuel elements allocated for 
heating the bleed gas. One axial and one radial section were used. 

Gas . - One axial section was used for both the energy and momentum equations of the 
gas in the bleed system. 


Nozzles 

Coolant . - Highly accurate nozzle-coolant calculations require many more than the 
one, two, or three sections feasible in analog work. Therefore, it was concluded that 
one section would be used to represent both the nozzle -coolant energy and momentum 
equations. During development of the nozzle model, exploratory analog calculations 
revealed that a one-section model yielded a momentum -pressure drop that was more 
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representative than either two- or three-section models. The comparison was made with 
results obtained from digital calculations in which 49 sections were used. 


* 

Turbopumps 

The turbine and pump characteristics were represented by appropriate functions set 
on variable diode function generators. Details are presented in appendix D. 


Valves 

The four bleed valves and the topping -turbine bypass valve were simulated by equa- 
tions presented in appendix D. Provision was made for both choked and nonchoked valve 
performance as applicable. 


Reactor Kinetics 

The neutron kinetics was represented by a one-energy group, space-independent 
model with six groups of delayed neutrons. Negative reactivity feedbacks were provided 
in response to the average temperatures of the fuel and core water. The positive contri- 
bution to reactivity of the core-gas density is negligible in the range of this investigation. 


VALIDATION 

After each computer was implemented with its part of the simulation, a static check 
was made to verify the scaling in the equations and computer patching. System condi- 
tions at the 100-percent design point were used in making the static checks. 

In addition to the 100-percent design point, conditions at two other arbitrarily se- 
lected reference points were determined for use in further checking of the model. These 
were at 60 and 36. 6 percent of design mass-flow rate and at design chamber temperature. 
Conditions at these three reference points (100, 60, and 36. 6 percent) were obtained 
from independent detailed digital and hand calculations made individually of the heat- 
exchanger, core, and turbopumps. As an aid in making the analog computations, con- 
trols were placed on poison reactivity and the topping -turbine bypass valve to set cham- 
ber temperature and core-gas mass-flow rate, respectively. Details of these controls 
will be presented subsequently. 
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TABLE II. - DEVIATION OF ANALOG RESULTS COMPARED WITH DIGITAL- AND 


HAND -CALCULATED RESULTS AT THREE REFERENCE POINTS 


Variables 

Percent of design value 

100 

60 

36.6 

Deviation, percent 

Reactor power 

- 0.93 

0 

0 

Hydrogen mass -flow rate through heat -exchanger, M x 

3 

0 

0 

Bleed-pump exit pressure, P pp E 

0.74 

0,84 

0 

Coolant -nozzle inlet pressure, j 

1 . 54 

1.6 

-14 

Coolant -nozzle exit pressure, P c ^ E 

0 . 55 

- 2 . 5 

- 2.6 

Hydrogen pressure at heat exchanger exit, P x E 

0 . 84 

- 1.1 

- 0.6 

Hydrogen temperature at heat exchanger exit, T x £ 

- 0.33 

0 

- 3.4 

Topping -turbine rotational speed, N TT 

- 0.11 

0 

- 11.8 

Topping -turbine mass -flow rate, fa TT 

3.1 

0 

- 14.5 

Mass -flow rate through topping -turbine bypass valve, M & 

- 11.9 

- 0 . 87 

15.1 

Hydrogen pressure at core inlet, P c j 

0 . 57 

- 0.37 

0.76 

Hydrogen temperature at core inlet, T c j 

- 0.31 

0 

- 3.3 

Hydrogen mass -flow rate in core, 

0 

0 

0 

Hydrogen pressure at core exit, P £ E 

- 0 . 5 

- 0 . 69 

- 0 . 22 

Hydrogen temperature at core exit, T c E 

0 

- 0.04 

0 

Bleed -gas mass -flow rate, M pG 

0 

0 

- 0.45 

Pressure at water -pump turbine throttle valve inlet, P 1 

- 1.2 

- 0 . 82 

- 3.4 

Temperature at water -pump turbine throttle valve inlet, T-^ 

1 

l -» 

OO 

- 0.43 

- 0.48 

Mass -flow rate through water -pump turbine throttle valve, M^ 

.16 

- 0.66 

- 1.7 

Mass -flow rate through water -pump turbine bypass valve, M 2 

- 0.63 

1.5 

1.2 

Water -pump turbine rotational speed, N WT 

0 

0 

0 

Pressure at water -pump turbine bypass valve exit, P 2 E 

- 2.0 

- 1.1 

- 3.7 

Temperature at water -pump turbine bypass valve exit, T 2 E 

- 1.7 

- 0 . 44 

- 0.49 

Mass -flow rate through bleed -pump turbine throttle valve, 

- 0.31 

- 0.23 

- 1.6 

Mass -flow rate through bleed -pump turbine bypass valve, M^ 

1.3 

- 0.18 

1.9 

Bleed -pump turbine rotational speed, 

- 0.18 

0 

0 

Pressure at bleed-pump turbine bypass valve exit, P^ E 

- 3.4 

- 5.6 

- 6.7 

Temperature at bleed -pump turbine bypass valve exit, E 

- 1.82 

- 0.47 

- 0.46 

Water mass -flow rate through heat exchanger, M x 

0 

0 

- 0.44 

Pressure at water-pump exit, P wp E 

0 

0 

0 

Mixed mean water temperature, T mmw 

0.06 

1 . 8 

- 0.23 

Water temperature at heat -exchanger exit, T x E 

- 0.28 

2.0 

- 0.05 



In order to check a reference point, the water-pump turbine throttle and bypass 
valves were first set at the reference positions. Then, minor adjustments (from the ref- 
erence settings) were made to the bleed-pump-turbine throttle and bypass valves until si- 
multaneous reference values were obtained for the core-gas exit temperature, core-gas- 
mass-flow rate, bleed-gas mass-flow rate, bleed-turbine mass-flow rate, and both 
bleed-turbine rotational speeds. Friction terms in the bleed system and water loop were 
trimmed slightly to attain the reference values. Steady-state data were recorded of all 
calculated parameters for comparison with the digital reference values. 

In table n (p. 11) are presented the percent deviations from the results of digital 
calculations of some of the more significant system variables at each of the three refer- 
ence points. It is immediately apparent that the only large deviations are in the topping 
turbopump, primarily at the 36. 6-percent point. In the initial setting of this point, it was 
observed that the turbopump variables agreed closely but the value of core -inlet pres- 
sure deviated slightly from the check point. Therefore, the friction in the lines joining 
the topping turbine with both the heat exchanger and core was adjusted to approach closely 
the reference value of core-inlet pressure. Friction settings that distributed the errors 
in these components evenly resulted in deviations of the order of 5 percent in each; how- 
ever, it was considered more important to the performance of the rest of the system to 
achieve reference conditions in the core. 

It is shown in table n that most of the other deviations throughout the system are of 
the order of 1 percent, with only two greater than 5 percent. A 5 -percent error is con- 
sidered acceptable for a nonlinear analog study of the magnitude undertaken herein. 


MAPPING PROCEDURE 

Three complete sets of steady- state maps were obtained. A set was made with the 
four bleed valves (water-pump-turbine throttle and bypass valves, and bleed-pump-turbine 
throttle and bypass valves) positioned at the 100-, 60-, and 36.6-percent reference points. 
These settings provided a topping-to-bleed-pump work split of 2 to 1 at each of the three 
reference points. The reason for maintaining fixed bleed- valves settings in the mapping 
survey was to determine if system control could be simplified by eliminating the need for 
these four manipulated variables. 

In order to expedite the data-taking process, proportional-plus-integral controls 
were placed on the poison reactivity and topping -turbine bypass valve, as was done in the 
validation study. The controls automatically regulated these two components to yield val- 
ues of core-gas exit temperature T „ (chamber temperature) and core-gas mass-flow 

. C, Uj 

rate M c in response to preselected demand signals. The controls were simulated by 
the following transfer functions: 
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For T 


c, E’ 


for M c , 


% 

e T 






, K 2,M 
s 


( 1 ) 


( 2 ) 


The error terms are the differences between demand and output signals; that is, 


and 


e T " ^ T demand ' T ^c, E 
e M = ^demand ' ^c 


(3) 

(4) 


With the mass-flow control set to maintain the topping -turbine bypass valve at a 
fixed position, a survey was made in the temperature range between 3000° and 6000° R 
by adjusting the demand signal of the chamber-temperature control to yield data at sev- 
eral selected values of T c g. After a traverse at constant hydrogen mass-flow rate, 
the topping -turbine bypass valve was repositioned by setting a different demand in the 
mass-flow control, and chamber temperature was again varied through its control. This 
procedure was repeated while the mass-flow rate was maintained constant at values of 
20, 30, 36.6, 40, 50, 60, 70, 80, 90, 100, 110, and 120 percent of the design value. The 
ranges of T £ g and M c were wide enough to provide data well beyond the design con- 
straints. In fact, as T c E was raised at a fixed mass-flow rate setting, four succes- 
sive constraints were exceeded. Thus, if one of a series of constraints should be re- 
laxed in the future, the margins from the next successive limitation would be available. 

The constraints applied to the mapping investigation of the system are as follows: 


Fuel temperature rise, °R 83 

Nozzle -wall temperature rise, °R 65 

Water-pump turbine-inlet temperature rise, °R 100 

Pressure -tube temperature rise, °R 20 

Heat -exchanger minimum water-side wall temperature, °R 492 

Bleed-pump stall, [(M/N)/(M/N) d 1 0. 723 

Turbine rotational speed rise, percent 10 
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RESULTS AND DISCUSSION 


The first important result of this study is that the system could be maintained at 
steady -state conditions with no external feedback controls. It was not known prior to 
operation of the simulation on the analog whether or not the engine system could be main- 
tained at steady-state conditions without closed-loop controls. 

The data obtained are plotted in figures 2 to 4 for bleed-valves settings correspond- 
ing to the 100-, 60-, and 36. 6 -percent reference points, respectively. All these figures 
are drawn over a background map of chamber pressure plotted against chamber tempera- 
ture with constant core-gas mass-flow rates. Chamber temperature is used as a coordi- 
nate because it is a primary factor in rocket specific impulse. Core-gas mass-flow rate 
and chamber pressure are representative of thrust. The variation in chamber tempera- 
ture on the maps is restricted to the range between 3000° and 6000° R to provide in- 
creased sensitivity. 

The background maps are all the same, irrespective of bleed-valves settings. Parts 

(a) present the locations of the various engine constraints. Parts (b) to (g) show curves 
of constant values of some of the more significant engine variables. In addition, con- 
straint boundaries pertaining to the appropriate bleed-valves settings are located in parts 

(b) to (g) of figures 2, 3, and 4. Four overlays are available on request for comparing 
the locations of the constraint boundaries obtained from one set of bleed-valves settings 
with those of the other two sets. A request card for the overlays is bound at the back of 
the report. Study of the equations in appendix D should prove helpful in understanding 
the behavior of the system and in following the trends brought out in the subsequent dis- 
cussions. 


Constraints 

100-percent bleed-valves settings . - Figure 2(a) reveals that with the bleed valves 
set at the 100-percent positions a feasible operating region extends from about 52 to 
120 percent of the design core-gas mass-flow rate with bleed-pump stall and speed limits 
forming the lower and upper boundaries. However, heat-exchanger freezing and the 
nozzle -wall temperature limit confine the allowable chamber temperature to a corridor 
between about 4000° and 4500° R. The margin of the design point from heat-exchanger 
freezing is greater than from the nozzle-wall temperature limit. 

It is evident from figures 2(a), 3(a), and 4(a) that the main hindrance to low-flow 
operation is the onset of stall in the hydrogen bleed pump. This observation indicates 
that as the total system mass-flow rate is reduced at constant bleed-valves settings, the 
bleed pump assumes a greater fraction of the load. 

Of primary significance in figure 2(a) is the intersection of the bleed-pump stall and 
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Core-gas exit temperature, T c £, °R 
(b) Reactor power. 

Figure 2. * Steady-state maps of tungsten water- moderated reactor preliminary 
reference system with 100- percent bleed-valves settings. 
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Figure 2. - Concluded. 

heat -exchanger freezing curves. This intersection, at about 52 percent of design core- 
gas mass-flow rate and 47 percent of design reactor power, establishes a lower limit to 
steady-state operation of the preliminary reference system with 100-percent bleed-valves 
settings. Thus, the attempt to simplify system control by fixing the four bleed valves at 
the 100-percent settings was not altogether successful. The simplification is feasible 
above approximately 47 percent of design power provided that the system is operated to 
maintain the chamber-temperature level within the 500° R range indicated in figure 2(a). 

Figure 2(a) reveals that neither the water-pump turbine -inlet temperature nor the 
pressure-tube temperature is limiting when the bleed valves are set at their 100-percent 
positions. Bleed-pump turbine rotational speed sets an upper bound to the operating re- 
gion, but at mass-flow rates considerably higher (120 percent) than the system 100- 
percent design point. The 110 -per cent bleed-pump rotational- speed constraint was 
specified to avoid exceeding the allowable stresses in the pump and turbine blades. 

60 -percent bleed- valves settings . - When a different set of bleed-valves positions is 
used, it is evident that the bleed system, water loop, and any parameters affected there- 
by will experience changes. Figure 3(a) presents the locations of the various constraints 
for the 60-percent bleed-valves settings. Comparison of figures 2(a) and 3(a) or placing 
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overlay 1 on figure 3(a) shows that changing the bleed-valves positions merely shifts the 
locations of the bleed-dependent constraints on the background map of chamber conditions. 
Only the nozzle-wall and fuel temperature limit curves are unaffected by the change in 
bleed-valves settings. Location of the nozzle-wall temperature limit curve is not altered 
because this temperature is assumed to be proportional to the chamber temperature (see 
appendix B and ref. 5). The fuel temperature is not influenced by the bleed system. For 
any given combination of core-gas exit pressure, temperature, and mass-flow rate, 
there is a corresponding value of fuel temperature. Any differences between the Tp 
curves of figures 2(a) and 3(a) are minor and result from slight changes in heat-transfer 
coefficients and gas properties. 

The most noticeable feature of the comparison of figures 2(a) and 3(a) or of overlay 1 
and figure 3(a) is the drastic reduction in the feasible operating region that results from 
adjusting the bleed valves to the 60-percent settings. The operating corridor is narrowed 
as the heat-exchanger icing curve is displaced to a position of higher chamber tempera- 
ture, while the water-pump turbine -inlet, and pressure-tube temperature limits are lo- 
cated at lower chamber-temperature positions. Displacement of the pressure-tube 
temperature-limit curve can be seen by comparing figures 2(a) and 3(a) or by placing 
overlay 2 on figure 2(a). 

On the other hand, closing down the bleed valves has the advantage of lowering the 
bleed-pump stall curve. At the 60-percent bleed- valves settings, the stall and heat- 
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exchanger freezing intersection occurs at about 39 percent of design flow and 37 percent 
of reactor power. As before, the bleed pump assumes a greater fraction of the load as 
the system mass-flow rate is lowered. When the bleed-pump-turbine throttle valve set- 
ting is reduced to the position corresponding to the 2-to-l pump power split at the 60- 
percent point, however, the load fraction taken by the bleed pump for any given set of 
chamber conditions is lowered. Therefore, the onset of bleed-pump stall is postponed 
during reduction of system mass -flow rate so that the stall line is lowered to the location 
shown in figure 3(a). 

The opposing effects of heat transfer to and from the water in the core and heat ex- 
changer (to be discussed subsequently) differ from those occurring at the 100-percent 
bleed-valves settings; thus, the heat-exchanger freezing curve is displaced to a more 
restrictive location in figure 3(a). 

At the reduced openings of the bleed valves, the bleed-flow rate is below that re- 
quired for adequate cooling of the bleed fuel assemblies at many points in figure 3(a) 
above a chamber temperature of 4200° R. The resulting rise in bleed-gas temperature 
causes the allowable water-pump turbine-inlet temperature to be exceeded, as shown by 
the latter curve cutting through the operating region of figure 3(a). Design chamber- 
temperature operation is unobtainable above approximately 67 percent of flow and power. 

Thus, the water-pump turbine -inlet temperature has replaced the nozzle -wall tem- 
perature as a primary constraint. An obvious way to avoid this restriction would be to 
open the water-pump-turbine and bleed-pump-turbine bypass valves to allow increased 
mass-flow rate through the bleed fuel elements. Doing so, however, would depart from 
the attempt to maintain constant bleed-valves settings. If, however, the limitation on 
T^ could be relaxed to permit turbine operation at 2050° R, most of the region between 
the heat-exchanger freezing and nozzle-wall temperature limits (including the 100- 
percent design point) would become available. Opening the operating region in this way 
to simplify control by fixing the bleed valves at the 60-percent settings would require 
advanced turbine technology, such as the use of higher-temperature-resistant materials 
or turbine cooling. 

Another interpretation of this observation is that in order for the bleed-pump tur- 
bine to develop the amount of work needed to achieve the 100-percent design point, its 
inlet temperature and therefore Tj must rise to 2050° R owing to the decrease in bleed- 
gas mass-flow rate. 

Closing down the bleed valves also decreases the water mass-flow rate at any point 

on the map in figure 3(a). Therefore, the pressure-tube temperature rises because of 

• 0 8 

the reduction in the heat -transfer coefficient caused by lowering the factor . in 
equation (D82). This explains the leftward displacement of the pressure-tube tempera- 
ture constraint on the map. 

36. 6-percent bleed-valves settings. - Examination of figure 4(a) reveals that an even 
smaller operating region is available when the 36. 6-percent bleed-valves settings are 
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used. The water-pump turbine-inlet temperature constraint cuts deeply through the cor- 
ridor bounded by heat-exchanger icing and the nozzle-wall temperature limit to elim- 
inate most of this envelope. The pressure-tube temperature limit also passes through 
this region. Meanwhile, however, the bleed-pump-stall curve is lower, and the heat- 
exchanger icing limit is relocated approximately at its position for the 100-percent 
settings. The nozzle -wall and fuel temperature limit boundaries are unchanged. 

With the 36. 6 -percent settings, the intersection of the bleed-pump stall and heat- 
exchanger icing limitations occurs at about 35 percent of design flow and 33 percent of 
power. In fact, the 36. 6-percent reference point itself lies on the bleed-pump stall line. 
Design chamber -temperature operation is restricted to mass-flow rates in the range be- 
tween 37 and 42 percent of design. The small operating region available with 36. 6- 
percent bleed- valves settings indicates that it would not be fruitful to operate with bleed- 
valves settings of less than the 60-percent positions. 

Comparison of figures 2(a) and 4(a) or placing overlay 1 on figure 4(a) shows that a 
larger operating region would be available with 36. 6 percent settings than with 100- 
percent settings if the water-pump turbine-inlet temperature limit could be relaxed. The 
lowered stall curve at 36. 6 percent accounts for this observation. It would be a distinct 
advantage in control if operation maneuvers throughout this region could be performed 
without requiring changes in the bleed-valves settings. To achieve the 100-percent de- 
sign point, however, a temperature of 2300° R is required at the water -pump turbine 
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inlet. Even if this were possible, 100-percent design-point operation would be denied by 
the pressure -tube temperature constraint. 


Composite Map 


In order to present the overall feasible operating region, a composite map of figures 
2(a), 3(a), and 4(a) was drawn. The resulting plot is shown in figure 5 as a plot of per- 
cent reactor power against hydrogen mass-flow rate. In the production of this figure, a 
continuous transition of constraint boundaries was assumed between figures 2(a) and 3(a) 
and between 3(a) and 4(a). Note that this plot is, in effect, merely a section of the bleed- 
pump head-flow performance map. Steady-state operation is confined to the enclosed re- 
gion, which is bounded by the bleed-pump stall curve at the lower end and by one of its 
speed curves at the upper end. A temperature limit in the nozzle wall bounds the upper 
side, and heat -exchanger freezing, the lower. 

The map in figure 5 summarizes the data presented in figures 2(a), 2(b), 3(a), and 
4(a). Its main advantage lies in providing an overall perspective of the capabilities and 
shortcomings of the preliminary reference design. The previously mentioned freezing- 
stall intersection establishing the lower bound of feasible steady-state operation at 33 per- 
cent of design reactor power is prominent. From this observation, the steady-state re- 
sults indicate that startup of the preliminary reference system may be a major problem. 

A conclusive determination of the startup capability of this system would require a de- 
tailed transient investigation in the low-power range. 



Figure 5. - Composite steady-state operating map of tungsten water- 
moderated reactor preliminary reference system. 


The margin of the 100 -percent de- 
sign point from heat -exchanger icing is 
much greater than that from the nozzle - 
wall temperature limit. Design 
chamber-temperature is represented, 
although not shown, by a 45° line 
through the origin in figure 5. This 
45° line lies just within the nozzle -wall 
temperature-limit boundary in figure 5. 
Thus, design specific impulse is poten- 
tially obtainable within the operating 
envelope above the 36. 6-percent refer- 
ence point. Although marginally ob- 
tainable in steady-state operation, de- 
sign chamber temperature may be un- 
obtainable from a practical standpoint 
where adequate margins are required. 
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The gradual narrowing of the operating corridor in figure 5 shows that in the preliminary 
‘reference system the heat-exchanger problem becomes more serious with decreasing 
system mass-flow rate. 


Engine Parameters 

Variations in some of the more important engine parameters are presented in fig- 
ures 2(b) to (f), 3(b) to (f), and 4(b) to (f). In these figures, constant values of each pa- 
rameter are plotted on the same background as that used in the discussion of constraints. 
Therefore, the superimposed constraint boundaries on the engine -parameter figures show 
the range of variation of the given parameter within the operating region. In the presen- 
tation that follows, most of the discussion is directed to figure 2 (100-percent bleed- 
valves settings) because it is typical of the values obtainable from figures 3 and 4. 

Reactor power. - The plot of figure 2(b) presents curves of constant reactor power. 
Comparison shows that figures 2(b), 3(b), and 4(b) are identical (except for the limit 
boundaries) despite the different bleed-valves settings. The value of reactor power is 
not influenced by bleed-valves settings, only by the combination of chamber conditions 
and core-gas mass-flow rate. These parameters correspond to a value of fuel tempera- 
ture, which in turn specifies 6kp, which together with 6k^ determines reactor flux and 
power. 

Within the operating envelope in figure 2(b), reactor power varies from about 47 to 
nearly 120 percent of design. Furthermore, each curve of constant power fraction 
crosses the corresponding mass-flow fraction line at the 100-percent design value of 
chamber temperature. This fact is explained readily by considering the following energy 
equation: 

Power = M c C P c (t e - T^ (5) 

At design core-gas exit temperature, or chamber temperature, the core-gas inlet tem- 
perature and specific heat are essentially constant. Therefore, reactor power and core- 
gas mass-flow rate assume the same fraction of their design values. 

Poison reactivity. - The relative locations of the operating envelope boundaries in 
figure 2(c) indicate the amount of poison reactivity required for the system to function in- 
side the operating region. The poison reactivity varies from 0 to about 1. 5 dollars within 
this enclosure. The 100-percent design-point value is about 35 cents, which does not in- 
clude any excess reactivity for xenon override. A variation in poison reactivity of about 
2 dollars characterizes the region between heat -exchanger freezing and the pressure-tube 
temperature limit. 

The curves of constant poison reactivity have a negative slope in figure 2(c). It 
would be desirable if they were vertical to simplify the controls. Then, if an operating 
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point was attained just above the pump stall line at a chamber temperature of 4460° R, 
the 100-percent design point could be reached by maintaining constant poison and merely * 
closing the topping-turbine bypass valve. A positive slope might even be better; a point 
in the lower left area of the operating region would be a starting point for ascent to the 
100-percent design point at constant poison reactivity. Nevertheless, knowledge of the 
location of the poison curves will be helpful in selecting an engine control system. Com- 
parison of figures 2(c), 3(c), and 4(c) reveals a displacement of the constant poison reac- 
tivity curves to the left at the lower bleed-valves settings, so that less poison is required 
at a given chamber condition. The corresponding slopes in the graphite -reactor system 
are positive. 

Water mass-flow rate . - The characteristics of the water mass-flow rate are pre- 
sented in figure 2(d). Values in the operating region vary from about 75 to 105 percent of 
the design value, a much smaller range than the reactor-power variation. Furthermore, 
the relative flatness of the water mass-flow rate curves in figure 2(d) indicates that this 
parameter is insensitive to variations in chamber temperature. 

Comparison of figures 2(d), 3(d), and 4(d) shows an upward displacement of the water 
mass-flow curves as the bleed-valves settings are lowered. This fact occurs because 
when the water-pump turbine throttle valve opening is reduced, the water-pump speed 
decreases; thus, the water mass -flow rate at any given point on the background map is 
lowered. 

Mixed-mean water temperature . - Figure 2(e) presents curves of mixed-mean water 
temperature at the reactor exit. Within the operating region, this parameter varies from 
620° to 720° R, approximately. The water system is pressurized to prevent boiling. 

Comparison of figures 2(b), (d), and (e) shows that, as reactor power rises, both 
water mass-flow rate and water-temperature level rise. The rise in power is accom- 
panied by rising water -pump turbine -inlet temperature Tj, which causes the rotational 
speed of this turbine to increase. Therefore, the water mass-flow rate rises. With 
rising reactor power, the bleed-pump turbine speed also increases, so that total hydro- 
gen mass -flow rate rises. 

In the heat exchanger, the transfer of heat from the hot water to the cold hydrogen is 
increased by the rises in the mass-flow rate of both fluids. Therefore, in the heat ex- 
changer there is a tendency toward decreased water temperature as a result of the in- 
creased fluid mass-flow rates caused by the power rise. 

Meanwhile, the rise in reactor power produces an increased heat input to the water 
in the core and the reflector. In addition, heat transfer from the pressure tube and re- 
flector material is enhanced by the increased water mass-flow rate. Thus, in the reac- 
tor, the rise in power has a tendency to increase the water-temperature level. 

Within the operating range, the water mass-flow rate is maintained at a relatively 
high level, varying from 75 to 105 percent of design (fig. 2(d)), while the reactor power 
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varies from 47 to 118 percent (fig. 2(b)). The high level maintained in the water mass- 
flow rate decreases the effects on the water in the heat exchanger and reactor. The two 
opposing tendencies in the heat exchanger and reactor account for the small change in 
water temperature experienced within the operating region. This explains why the water 
temperature rises only about 100° Rwhen the reactor power level is raised from 47 to 
118 percent of design. 

Percent bleed . - The variation in percent bleed is presented in figure 2(f). For 100- 
percent design-point operation, almost 2. 6 -percent bleed is required. Within the oper- 
ating region, this parameter varies from less than 2. 5 to approximately 2. 8 percent. 

The percent bleed rises with decreasing total hydrogen mass-flow rate because the con- 
stant bleed-valves settings force the bleed system to pass more hydrogen gas than is 
needed at low system mass -flow rates. 

This trend of rising percent bleed with decreasing system flow rate also acts to 
maintain the variation in water mass-flow rate and temperature level in the small range 
previously noted. As reactor power is decreased, the greater-than-required percent 
bleed causes the water-pump turbine to rotate faster than required. Therefore, the water 
mass-flow rate is maintained at a higher value, and, although the water -temperature 
level is lowered because of the power decrease, the drop is small because of the favor- 
able heat -transfer conditions in the reactor. 

A similar situation occurs when the corresponding plots for the 60- and 36. 6-percent 
bleed-valves settings are considered. Figures 3(f) and 4(f) show that at the 100-percent 
design point approximately 2. 3 and 2. 0 percent bleed occur, respectively. In these in- 
stances, the bleed-valves openings are closed more than required for high-power opera- 
tion. The decrease in percent bleed at any point (such as the 100-percent design point) 
results in higher bleed-gas temperature. This explains why the curve of limiting water- 
pump turbine-inlet temperature moves into the operating corridors (see figs. 2(a), 3(a), 
and 4(a)) as the bleed-valves settings are changed progressively from the 100 to 60 to 
36. 6 percent settings. 

Because of choking in the bleed valves, it is not necessary to control the bleed-gas 
mass-flow rate. With constant bleed-valves settings, the variation in percent bleed is 
small. A feature of the bleed system is its self -regulating characteristic: when hydro- 
gen mass-flow rate is decreased, the percent bleed increases. This fact is also indi- 
cated by the negative slope of the water-pump turbine -inlet temperature constraint in 
figures 2(a), 3(a), and 4(a). This adds to the safety of bleed fuel-element operation. 
There is a slight penalty incurred, however, in the reduction in turbine energy available 
at low hydrogen mass-flow rates. 

Topping-turbine bypass valve setting . - During the mapping survey, the position of 
the topping-turbine bypass valve was recorded. From this information, lines of constant 
valve setting are presented in figure 2(g), with the value at the 100-percent design point 
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as a base. These lines are similar to the core-gas mass-flow rate lines, except that 
their slopes are steeper. Their dissimilarity is due to two main factors: friction in 
core gas and topping -turbine inefficiency. The curves of constant valve setting and 
mass-flow rate are most dissimilar at the high-flow region in figure 2(g). In this re- 
gion, the topping-turbine efficiency is close to its design value; thus, this factor has a 
small effect. However, the friction pressure drop, which is proportional to M^, is 
high so that the valve opening is increased over that required for frictionless flow. In 
the lower regions of figure 2(g), a low turbine efficiency resulting from far-off -design 
turbine operation is a major factor. Frictional pressure drop has a smaller effect 
here. 

Along a line of constant M c in figure 2(g), increasing chamber temperature re- 
quires opening the topping -turbine bypass valve, and closing the valve is accompanied by 
decreasing temperature. This means that at constant core-gas mass-flow rate the top- 
ping turbine develops a smaller share of the total work as chamber temperature is in- 
creased. Accordingly, the bleed-pump turbine produces a larger share of the work. 


OPERATION AND CONTROL 

In addition to defining the boundaries of feasible system operation, an objective of 
this study was to explore the system operation and control characteristics. The mapping 
surveys were conducted with the use of fixed bleed-valves settings to determine if this 
simplification in control could be used over a wide operating range. Figure 2(b) shows 
that this is possible above approximately 47 percent of reactor power with the use of the 
100-percent bleed-valves settings. Of the six manipulated variables for control listed on 
page 6, this simplification requires that only two, poison reactivity and topping-turbine 
bypass valve, be adjusted for maneuvers within the operating envelope of the map given 
in figure 2(a). 

On the other hand, the severely restrictive confines of the operating corridor in fig- 
ure 5 indicate that simplification in control may not be a valid objective. Because ma- 
neuvers within this envelope are so limited, it may be advisable to monitor the various 
damage limits and provide the capability to override the main control parameters. In 
doing so, however, measuring such parameters as nozzle-wall temperature accurately 
may introduce a major problem. 

Since these practical considerations have been recognized, this report discusses the 
theoretical aspects in further simplifying system control. The ultimate simplification 
would arise if the need for controlling either poison reactivity or the topping -turbine by- 
pass valve could be eliminated. Thus, it would be highly desirable to control system 
operation by the topping -turbine bypass valve alone. The nonfeasibility of this procedure 
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is readily displayed by considering figure 2(c). From the 100-percent design point if 
constant poison reactivity is maintained, a maneuver involving opening the topping- 
turbine bypass valve results shortly in overheating the nozzle wall. The reverse proce- 
dure started from a lower point steers the system toward heat -exchanger icing. The 
negative slopes of the constant poison reactivity curves prevent this simplification. 

On the other hand, figure 2(g) reveals that fixing the topping -turbine bypass valve 
setting to attempt operating maneuvers on the map by adjusting poison reactivity alone is 
also in vain. Thus, with only one control, it is impossible to operate the system within 
the feasible operating region. A minimum of two control variables is required. 

Possible Operating Paths 

This discussion implies that two control parameters are required for maneuvers 
within the feasible operating envelope. Possible use of parameters other than poison 
reactivity and topping-turbine bypass valve is discussed in the next section. For the 
present, however, a simple path employing both of these control parameters is consid- 
ered. Possible modifications to the system that would permit startup and acceleration 
through the low -power region to the 100-percent design point are disscussed in the sec- 
tion SUGGESTED MODIFICATIONS. It is suggested that this method be used to attain 
a point within the operating capability of the preliminary reference system with the 



Figure 6. - Suggested operating path for tungsten water -moderated reactor system. 
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100-percent valves settings. Such a point is designated point A in figure 6 and on over- 
lay 4, which will be helpful in illustrating a path from A to the 100-percent design point D: 
Figure 6 consists of a background map of hydrogen mass-flow rate curves, a superim- 
posed operating region for 100 -percent bleed-valves settings, and a suggested operating • 
path taken from figures 2(c) and (g). Overlay 4 can be used with these two figures. At 
point A, the heat -exchanger splitter valve of the modified reference system can be closed, 
and regulation of the bleed-pump-turbine throttle and bypass valves can be discontinued 
by setting them at the 100-percent positions. The modified reference system therefore 
reverts to the preliminary reference system. 

Point A is characterized by the following system parameters: 

Powe i: = 0.56 

Design power 


T c , E = 41 8°° R 

8kl = $1.35 



If this topping-turbine bypass valve setting is maintained constant while the poison reac- 
tivity is reduced to 1 dollar, point B is reached. Thereupon, this valve can be closed 
down to the 100-percent design-point setting while the poison reactivity is kept at 1 dollar. 
This establishes point C. The 100-percent design point D is now accessible by decreasing 
the poison while keeping the topping-turbine bypass valve setting fixed at the 100-percent 
design position. The path ABCD is also included on the composite map (fig. 5, p. 28) so 
that the changes in reactor power and hydrogen mass-flow rate can be shown. From 
B to C, the reactor power rises from 60 to 90 percent of design, while the mass-flow 
rate rises from 60 to 97 percent. 

Thus, for path ABCD only two controls are needed, poison reactivity and topping- 
turbine bypass valve, only one of which is operated at a time. Of course, a more direct 
path from A to D is possible by controlling both of these parameters simultaneously. 
Furthermore, many alternate paths are possible within the limited operating envelope. 

If the constant poison reactivity curves had a positive slope, it might be possible to pro- 
ceed directly from A to D. Further study is required to determine if changes to the 
system design could yield a vertical or positive slope. 

From point B to C, the chamber temperature decreased from a value of 4450° to 
4120° R. If this drop should be judged unacceptable from the standpoints of vehicle per- 
formance or fuel-temperature considerations, constant -temperature paths such as A to 
C or B to D could be programed. 
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Possible Control Parameters 
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Figure 7. - Variation in reactor power with hydrogen mass- 
flow rate at constant poison reactivity of 50 cents. 
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Figure 8. - Variation in reactor power with mixed-mean water 
temperature. 


Hydrogen mass -flow rate . - In the 
search for a parameter suitable for a re- 
actor power controller, the hydrogen 
mass -flow rate was considered. For this 
purpose reactor power was plotted in fig- 
ure 7 against hydrogen mass-flow rate, 
with poison reactivity maintained constant 
at 50 cents for the three bleed-valves set- 
tings. Use of a 50-cent poison reactivity 
in figure 7 accounts for the fact that 
slightly less than 100-percent design re- 
actor power is indicated at 100-percent of 
design mass -flow rate. As was mentioned 
previously, the poison reactivity at the 
100-percent design point is about 35 cents 
(see fig. 2(c), p. 16). Figure 7 shows that 
reactor power is nearly proportional to 
hydrogen mass-flow rate. This relation 
is also true if chamber temperature is 
kept constant instead of poison reactivity 
(fig. 2(b), p. 15). Therefore, the topping- 
turbine bypass valve, through its effect on 
the hydrogen mass-flow rate, appears 
promising as a parameter to control reac- 
tor power. Because engine thrust is also 
proportional to hydrogen mass-flow rate, 
its effect on engine thrust must also be 
considered. 

Water temperature . - When this re- 
actor system was conceived, it was anti- 
cipated that water temperature might be 
used to control reactor power. For ex- 
ample, the low temperature of the water 
could be sensed more easily than hot-gas 
temperature. Figure 8 was plotted to ex- 
plore this possibility. In figure 8(a), 
relative reactor power is plotted against 
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mixed-mean water temperature at design chamber temperature for all three bleed-valves 
settings. The steep slopes of the curves in this figure reflect the previously mentioned • 
small change in water temperature accompanying large changes in reactor power. Fig- 
ure 8(a) indicates that, when the reactor power changes 1 percent of the design value, the 
water temperature changes about 1° R for the 100-percent bleed- valves settings. For the 
other two bleed-valves settings, the power change is about 0. 7 percent per °R. 

Figure 8(b) is similar to figure 8(a) except that poison reactivity is maintained con- 
stant at a value of 50 cents, and, therefore, chamber temperature varies. Also plotted 
in figure 8(b) are curves of constant fuel temperature change to aid in displaying the feed- 
back reactivity characteristics of the system. Under the condition of constant poison re- 
activity, less than a 3-percent change in reactor power is accompanied by a 1° R change 
in water temperature. Therefore, controlling reactor power by sensing water tempera- 
ture during either constant temperature or constant poison operation appears feasible. 

Figure 8(b) shows that along any of the curves of constant bleed-valves settings the 
water temperature and reactor power both rise (or drop) together; that is, the slopes of 
the curves are positive. That this is not a paradox, as it might seem from the negative 
water-temperature reactivity of the system, is shown by considering the relative loca- 
tions of the 'constant fuel-temperature change curves. As reactor power is increased 
along a curve of constant bleed-valves settings in figure 8(b), the fuel-temperature level 
decreases. This decrease in fuel temperature causes an increase in the fuel reactivity 
5kp (see eq. (D178)), which tends to raise reactor power. At the same time, the de- 
crease in 6k^ in response to the increase in water temperature tends to lower the reac- 
tor power level. Under the condition of constant poison reactivity in figure 8(b), the re- 
sulting reactivity is determined by 6k' F and 6k^ (see eq. (D176)) only. The positive 
slopes shown in figure 8(b) imply that the effect of 6kp is dominant under these condi- 
tions. This observation follows from the fact that, when the fuel temperature drops 
200° R, in figure 8(b), the water temperature rises only 6° or 7° R. 

The constant fuel-temperature change curves in figure 8(b) provide further insight in 
understanding the behavior of the system with respect to temperature-reactivity feed- 
backs. Because poison reactivity is maintained constant throughout figure 8(b), the 
curves of constant fuel-temperature change isolate the effect of the water -temperature 
reactivity on reactor power. Equation (D178) shows that, if fuel temperature is constant, 
6kp is constant. Therefore, in accordance with equation (D176), 5k^ alone is respon- 
sible for the reactor power variation along the curves of constant fuel-temperature 
change. Thus, from the curve representing 100-percent bleed-valves settings if the 
system is operated to maintain constant fuel temperature, figure 8(b) discloses that the 
water temperature rises. The concomitant decrease in 5k^ results in a drop in reactor 
power. This hypothetical manner of operating the system necessitates changes in the 
bleed-valves settings. 
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The negative slopes of the constant fuel- 
temperature change curves reflect the nega- 
tive characteristic of the water -temperature 
reactivity. The steepness of the slopes illus- 
trates the high sensitivity of reactor power to 
this parameter. The safety aspect inherent in 
this reactor design concept is shown clearly 
by the fact that the slopes of the constant fuel- 
temperature change curves become steeper 
with increasing water temperature. 

Water mass -flow rate. - Figure 9 is a 
cross plot made to present the effect of water 
mass-flow rate on reactor power when poison 
reactivity is maintained at 50 cents. Once 
again, curves of constant fuel-temperature 
change are included to enhance the understanding of interactions caused by temperature- 
reactivity feedbacks. At the 100-percent bleed-valves settings, a 2-percent change in 
reactor power is accompanied by a 1-percent change in water mass-flow rate. The sen- 
sitivity of water mass-flow rate to reactor power variation implies that sensing this flow 
rate might be employed for reactor-power control. Other factors, however, outweigh 
this characteristic. Water mass -flow rate should be maintained at a high level to allevi- 
ate the possibility of heat -exchanger icing and to provide an extra margin from the 
pressure-tube temperature limit. Furthermore, a high water mass-flow rate is indica- 
tive of a fast water-loop response time, which is desirable because it increases the tem- 
perature coupling of the hydrogen and the water. 

Figure 9 shows that along any curve representing constant bleed-valves settings and 
with constant poison reactivity, water mass-flow rate rises when reactor power rises. 
The water mass-flow rate rises despite a drop in percent bleed (fig. 2(f)). This drop is 
small so that the large increase in total hydrogen mass-flow rate (fig. 2(c)) results in a 
rise in bleed-gas mass-flow rate and, hence, a rise in water-pump turbine energy. The 
increase in water mass-flow rate results in raising the water-temperature level but only 
by a small amount, as explained in the discussion of the opposing effects in the heat- 
exchanger and reactor. Finally, as discussed in connection with figure 8(b), the effect 
of the small water-temperature rise is outweighed by the 200° R drop in fuel-temperature 
so that a net rise in reactor power is experienced. 

As in figure 8(b), the constant fuel-temperature change curves in figure 9 isolate the 
effect of the water. In this instance, however, the effect is interpreted in terms of water 
mass-flow rate because of its use as the abscissa in figure 9. The positive slopes of the 
curves of constant fuel-temperature change (fig. 9) show that reducing the bleed-valves 



Figure 9. - Variation in reactor power with water mass- 
flow rate at constant poison reactivity of 50 cents. 
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settings results in decreased water mass-flow rate and reactor power. In the same hypo- 
thetical mode of operation, figure 8(b) indicates that the water temperature rises. This . 
rise accounts for the observation in figure 9 that reactor power decreases along a con- 
stant fuel-temperature change curve as the bleed-valves settings are reduced. 

Poison reactivity . - Comparison of figures 2(c) and (e) reveals that the slopes of the 
constant poison reactivity and water-temperature curves at any given point are similar. 
The slopes of the former curves are somewhat steeper, however. The similarity is not 
surprising when it is recalled that constant fuel-temperature curves are nearly vertical 
in figure 2(a). This characteristic lessens the comparative effect of Sk'p, on the curves 
of figures 2(c) and (e). Thus, because the net reactivity is zero in steady-state operation, 
equation (D176) implies that 6k^ and 6k^ tend to offset each other. That is, the con- 
stant poison- reactivity curves are a direct reflection of the water-temperature reactivity, 
which, in turn, follows the water temperature. 

The results implied in the foregoing discussion are verified by figure 10, which 
shows the variation in reactor power with poison reactivity when chamber temperature is 
maintained at 4460° R. Because the poison reactivity is essentially the negative of the 
water reactivity at constant chamber temperature, figure 10 is approximately a mirror 
image of figure 8(a). The lack of exact correspondence lies in the fact that the fuel reac- 
tivity varies slightly along each of the curves. As in figure 8(a), the steepest slope oc- 
curs with 100-percent bleed-valves settings. At the latter settings, insertion of 10 cents 
of poison results in a decrease of about 7 percent in the design power. Furthermore, 
with poison insertion there is a transient tendency to undershoot both reactor power and 
core -gas exit temperature. 



Figure 10. - Variation in reactor 
power with poison reactivity at 
design chamber temperature of 
4460° R. 


The extreme sensitivity of reactor power to poison re- 
activity requires that reactor power be accurately sensed 
and that poison reactivity be minutely controlled for opera- 
tional maneuvers, such as path ABCD in figure 5. In the 
modifications to be discussed subsequently, it was found that 
the heat -exchanger splitter valve was suitable for control of 
reactor power. If it should be found that finer control is 
available from manipulation of this valve, once point A in 
figure 5 is reached by the modified reference system, to 
reach point D it may be desirable to make the appropriate 
changes in reactor power by adjusting the heat-exchanger 
splitter valve instead of the poison reactivity. Furthermore, 
the tendency of poison insertion to cause transient tempera- 
ture undershoots would be replaced by a smoother transient 
resulting from the transport delay in the water loop. Per- 
formance obtained by adjusting the heat -exchanger splitter 
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valve with constant poison reactivity would approximate that using water temperature as 
a controller. 


Bleed-Valves Considerations 

It has been shown that the threshold of bleed-pump stall recedes to a lower hydrogen 
mass-flow rate as the bleed-pump-turbine throttle valve setting is reduced to satisfy the 
requirements of a 2-to-l power split between the topping and bleed pumps. This obser- 
vation indicates that fixing the bleed-pump-turbine throttle valve at a minimum setting 
would yield the lowest stall condition. Then, since the system tends to unload the bleed 
pump as mass-flow rate and power rise, the bleed-pump power fraction would be at a 
minimum at the 100-percent point. Because of the nature of pumps in series, the topping 
pump assumes the load fraction dropped by the bleed pump and thus causes the former to 
approach stall. There evidently exists a problem of optimizing the fixed bleed-pump- 
turbine throttle valve setting to accommodate the stall limits of both pumps. It does not 
appear necessary to control the pump work split in the steady-state range above 47 per- 
cent of design reactor power because of the absence of stall above this level. 

Because of the trend for the bleed pump to increase its load fraction as flow and 
power decrease, it appears that this pump would carry the full load at startup. Refer- 
ence 5 corroborates this expectation which, however, cannot materialize because of the 
pump stall characteristic. Thus the topping pump must carry a share of the load. There- 
fore, a control would be required to maintain a suitable, varying power split during 
startup. An alternate scheme would be to provide a control to maintain operation of both 
pumps with equal margins from stall, as was done in references 5 and 6. 

A high water mass-flow rate is desirable to avoid heat -exchanger freezing and for 
decreasing the water-loop response time. Therefore, it is concluded that both the water- 
pump-turbine throttle and bypass valves can be replaced by orifices sized at the 100- 
percent conditions for use over the operating range of this investigation. Sizing the ori- 
fices at the 100-percent conditions would provide adequate water mass-flow rate and 
would prevent the water-pump-turbine-inlet temperature limit being reached at low re- 
actor power levels. 

The bleed-pump-turbine throttle and bypass valves cannot be replaced by orifices 
for operation below about 47 percent of design reactor power if they are set to yield a 
2-to-l pump work split at the 100-percent design point because of bleed-pump stall. With 
the heat-exchanger bypass configuration, these two valves require continuous manipula- 
tion to permit startup and acceleration through the low -power range. 


SUGGESTED MODIFICATIONS 


It is evident from figure 5 (p. 28) that improvement in performance of the prelimi- 
nary reference system requires that the operating envelope be opened to encompass a 
wider range and extended to allow low-power operation. The location of the high-power 
side of the operating corridor is considered critical because although it does not theoreti- 
cally restrict design-chamber-temperature operation, it may not provide an adequate 
margin from the nozzle -wall temperature limit. Movement of this and the successive 
boundaries would require relaxing the nozzle wall, fuel, and water-pump turbine-inlet 
temperature limits or changing the design of the nozzle and the main bleed fuel elements. 

Widening the operating region also involves finding a means of overcoming the heat- 
exchanger freezing problem to allow removal or displacement of this boundary on the 
map of figure 5. To do so requires reducing the water -temperature drop in the heat ex- 
changer. The drop in water temperature can be reduced by increasing the water mass- 
flow rate, decreasing the exchanger surface area[, or decreasing the hydrogen mass- 
flow rate through the heat exchanger. Raising the water mass-flow rate would be ineffec- 
tive for this purpose, however, because of the opposing effects in the reactor and heat 
exchanger on water-temperature level. Reducing the heat-exchanger surface area enough 
to be effective would cause the pressure-tube temperature limit to be exceeded because 
of inadequate cooling. Decreasing the hydrogen mass -flow rate in the heat exchanger 
also tends to promote excessive pressure-tube temperature because of the resultant 
higher heat -exchanger water-exit temperature. This scheme, however, has the advan- 
tage of being controllable. A portion of the main hydrogen mass flow can be bypassed 
around the heat exchanger, the amount being regulated by a valve. 

This discussion implies that heat -ex changer icing and a pressure-tube temperature 
limit are fundamental limits in the water system. Any modification selected must pro- 
vide a compromise that constrains the water system to operate between these two limits. 

Extending the operating envelope farther into the low-power region requires lessen- 
ing the load on the bleed pump to keep it out of stall. Programing control of the bleed 
valves to maintain both bleed and topping pumps equidistant from stall is a possibility. 
Decreasing the hydrogen system pressure drop, or impedance, is another. This could 
be accomplished by bypassing some of the hydrogen around the nozzle-coolant tubes, 
which form a high-pressure drop component. 

In accordance with these considerations, three particular modifications to the pre- 
liminary reference system were selected for further study. This work is reported in 
references 5 and 6. 

In one modification, the preliminary reference system was altered to include hydro- 
gen bypass around the heat -exchanger with a splitter valve at the exit. It was shown that 
this modification successfully enlarged the operating region and extended low-power 
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operation down to about 13 percent of design reactor power. A second modification in- 
cluded nozzle -coolant tube bypass in addition to the he at -ex changer bypass. This modifi- 
cation was inferior to the preliminary reference system in high-power performance. 
Emergence of the nozzle-coolant temperature-limit constraint prevented design chamber- 
temperature operation over the entire range below the 100-percent design point. The 
third modification studied was an all -topping cycle with heat -exchanger hydrogen bypass. 
Besides achieving virtually the same performance improvement as the first modification, 
this design also simplified controls and system complexity by elimination of the bleed 
system. 

From the results of references 5 and 6, it appears that system startup would be pos- 
sible by using any of these modifications. The all-topping cycle is ideal from the stand- 
point of requiring control of only the topping -turbine and heat-exchanger bypass valves. 
The topping cycle, however, appears to have limited growth potential for developing the 
high chamber pressures anticipated for use with advanced nozzle designs. This limitation 
arises from the inability of the topping cycle to supply the turbine with enough energy to 
meet increased demands in turbine power within the existing core design. The modifica- 
tion using heat -ex changer bypass also has startup capability, but requires continuous 
regulation of the bleed-pump-turbine throttle and bypass valves to maintain both hydrogen 
pumps at equal margins from stall. This modification does, however, allow replacement 
of the water -pump-turbine throttle and bypass valves by fixed orifices. The heat- 
exchanger bypass modification was therefore chosen for more extensive study. With re- 
gard to the preliminary reference system again, it may be possible to preheat the water 
to permit chilldown and startup without heat -exchanger icing. 

As a consequence of results obtained from the study of the preliminary reference 
system, a simple modification to the design incorporating a cold-hydrogen bypass around 
the heat-exchanger was investigated in references 5 and 6. Some of the conclusions made 
by references 5 and 6 pertaining to the modified reference system are as follows: 

(1) The operating corridor was opened considerably, and a limit on the pressure- 
tube temperature replaced heat-exchanger freezing as the low chamber-temperature 
boundary. 

(2) Startup appears feasible, because the steady-state pump-stall-pressure-tube 
temperature-limit intersection was lowered to about 13 percent of design reactor power. 

(3) Control of the heat -exchanger splitter valve can function as an alternate or 
supplement to poison reactivity. 
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CONCLUSIONS 


The steady -state analog simulation of the tungsten water -moderated reactor prelim- 
inary reference system led to many conclusions. Some of the more significant conclu- 
sions are as follows: 

1. The system appears stable without external feedback controls. 

2. A lower limit to steady-state operation is imposed by the simultaneous occurrence 
of heat -exchanger icing and bleed-pump stall. This occurs at about 35 percent of design 
hydrogen mass-flow rate and 33 percent of design reactor power. 

3. Although the steady-state results indicate that the preliminary reference system 
may encounter severe problems in startup, conclusive evidence would require a low- 
power transient investigation. 

4. High-power operation is limited by bleed-pump turbine rotational speed. 

5. As chamber temperature is increased, the system encounters successive temper- 
ature limits in the following components: nozzle-wall, fuel, water-pump turbine -inlet, 
and pressure -tube. 

6. Reactor power can be controlled by two parameters: the topping -turbine bypass 
valve and poison reactivity. As an alternative, water temperature is a suitable control 
parameter. 

7. In a search for parameters suitable for reactor control, the following relations 
were found: reactor power is approximately proportional to hydrogen flow rate; a poison 
reactivity insertion of 10 cents reduces reactor power by about 7 percent; and less than 
a 3-percent change in reactor power is accompanied by a 1° R change in water- 
temperature level. 

8. The water mass -flow rate should be maintained at a high level; therefore, the 
water-pump-turbine throttle and bypass valves can be replaced by fixed orifices sized 
for 100 percent design-point operation. 

9. Reactor power and hydrogen mass-flow rate are not influenced by pump work 
split or the positions of the valves in the bleed system. 

10. In the operating range above 47 percent of design reactor power, system control 
can be simplified by maintaining the bleed-pump-turbine valves at their 100-percent 
positions. In this region, operation can be controlled by the topping -turbine bypass 
valve and poison reactivity, only. 

11. As reactor power is reduced with constant bleed- valves settings, the bleed pump 
assumes an increased fraction of the load causing this pump to approach stall. Below 
47 percent of design reactor power, the pump work split must be adjusted to keep the 
bleed pump out of stall. 
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12. The bleed-gas mass-flow rate does not require control in the range of this study. 
There is a self -regulating tendency in which the percent bleed rises as reactor power and 
total hydrogen mass flow rate are decreased. 

13. Design chamber temperature (or specific impulse) is potentially obtainable at all 
reactor power levels above about 37 percent of design, at which point bleed-pump stall 
occurs. The margin from the nozzle-wall temperature limit, however, is very small. 

14. The settings of the valves in the bleed system determine the locations of all sys- 
tem constraints on a steady-state map, except the nozzle -wall and fuel-temperature 
limits. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 15, 1966, 

122-28-02-04-22. 
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APPENDIX A 


SYMBOLS 


A 

area, ft^ 

BTBV 

bleed-pump-turbine bypass 
valve 

BTTV 

bleed-pump-turbine throttle 
valve 

b 

hydrogen temperature -ratio 
factor in expression for 
convection heat -transfer 
coefficient 

C i’ 

concentration of delayed- 
neutron precursors of i’* 11 

O 

group, nuclei/ft 

C P 

specific heat at constant pres- 
sure, Btu/(lb mass)(°R) 

Cy 

specific heat at constant vol- 
ume, Btu/ (lb mass)(°R) 

c v 

valve flow -capacity coefficient 

D 

hydraulic diameter, ft 

E 

thermodynamic internal en- 
ergy, Btu/lb mass 

e 

total internal energy, 
Btu/lb mass 

F' 

mass-flow function set on 
variable diode function 
generator 

f 

friction factor 

f' 

pressure -ratio function set 
on variable diode function 
generator 


f" water-temperature reactivity 

function set on variable 
diode function generator 

G factor in expression for heat- 

transfer coefficient 

G* pump pressure-speed function 

set on variable diode func- 
tion generator 

g acceleration of gravity, 

ft/sec^ 

g' pump torque -speed function 

set on variable diode func- 
tion generator 

H factor in expression for con- 

vection heat -transfer coef- 
ficient 

h convection heat -transfer coef- 

ficient, Btu/(ft^)(sec)(°R) 

I turbopump moment of inertia, 

lb mass (ft)(sec^) 

J mechanical equivalent of heat, 

(ft)(lb mass)/Btu 

K, Kj, Kg constants 

k thermal conductivity, Btu/ 

(ft)(sec)(°R) 

k' reactivity 

L torque, (lb force) (ft) 

l * mean lifetime of neutrons, 

sec 
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M 

mass-flow rate, lb mass/sec 

A, 

fraction of fission neutrons that 

N 

rotational speed, rpm 

l 

are delayed neutrons in i'** 1 

n 

P 

3 

neutron density, neutrons/ft 

o 

absolute pressure, lb force/ft 

r 

group 

ratio of specific heats for hydro- 
gen gas, Cp/C v 

Q 

heat generation rate, Btu/ sec 

€ 

error between demand and calcu- 

R 

gas constant for hydrogen, 


lated values 


(ft)(lb force)/(lb mass)(°R) 

e 

topping-turbine bypass valve 

Re 

Reynolds number 


angle, deg 

S 

general heat -flow rate, including 

V 

decay constant of delayed-neutron 

s 

generation and convection, 
Btu/ sec 

Laplace transform variable 

i 

M 

precursors of i' 1 group, 
sec - * 

absolute viscosity, lb mass/ 

T 

absolute temperature, °R 


ft- sec 

TTBV 

topping -turbine bypass valve 

P 

mass density, lb mass/ff 3 

t 

time, sec 

Subscripts: 

u 

fluid velocity, ft/ sec 

A 

upstream axial section in three - 

V 

WTBV 

volume, ft^ 

water -pump -turbine bypass 

B 

section model 

middle axial section in three- 
section model 

WTTV 

valve 

water-pump-turbine throttle 
valve 

B, E 
BF 

bleed exit 
bleed fuel 

X 

axial distance, ft 

BG 

bleed gas 

y 

radial distance, ft 

BN 

bleed nozzle 

z 

constant entrance-exit loss fac- 

BP 

bleed pump 


tor in expression for friction 

BT 

bleed-pump turbine 


factor 

BTP 

bleed turbopump 

z 

distance in plane normal to both 

b 

bulk 

0 

x and y, ft 

sum, for all groups, of frac- 

C 

downstream axial section in 
three- section model 


tions of fission neutrons 
that are delayed neutrons 

CN 

coolant nozzle 
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CNW 

coolant-nozzle wall 

S 

surface 

c 

hydrogen in core 

SH 

stagnant hydrogen 

D 

flow -divider tube 

T 

temperature 

des 

design 

TOT 

total 

E 

exit 

TP 

topping pump 

F 

fuel 

TT 

topping turbine 

f 

film 

TTP 

topping turbopump 

fr 

friction 

W 

water- side surface in heat ex- 

H 

hydrogen-side surface in heat 


changer 


exchanger 

WP 

water pump 

He 

heat exchanger 

WT 

water-pump turbine 

I 

inlet 

WTP 

water turbopump 

i 

inner 

w 

water in reactor 

i' 

index representing energy group 

X 

hydrogen in heat exchanger 


of delayed neutrons 

X 

water in heat exchanger 

M 

mass-flow rate 

1 

water-pump-turbine throttle valve 

max 

maximum 

2 

water-pump-turbine bypass valve 

mmw 

mixed mean water 

3 

bleed-pump-turbine throttle valve 

N 

nozzle hot-gas side 

4 

bleed-pump-turbine bypass valve 

0 

outer 

5 

topping-turbine bypass valve 

P 

pump 

1’ 

inlet axial station 

PT 

pressure tube 

2’ 

axial station at interface between 

P 

poison 


sections A and B 

R 

side reflector 

3’ 

axial station at interface between 

ref 

reference 


sections B and C 
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APPENDIX B 


ASSUMPTIONS 

The following assumptions were used in the analysis: 

General: 

(1) Flow through piping is adiabatic. 

(2) Temperature transport lags in piping are negligible. 

(3) Radiation heat transfer to the pressure tubes from the tungsten support tubes is 
included in the term for conduction through the stagnant hydrogen. 

(4) Gamma and neutron heating rates are proportional to reactor power level. 

(5) Radial temperature gradients in the fuel, pressure tubes, and flow dividers are 
negligible. 

(6) Heat loss entering the pressure tubes through the stagnant hydrogen is a linear 
function of the reactor power. 

(7) Friction factors are constants, but a different constant is used with each set of 
bleed- valves settings. 

Hydrogen Flow System: 

(1) Flow through both hydrogen pumps is liquid, and density is constant in each pump. 

(2) Hydrogen gas is treated as a perfect gas. 

(3) Gas properties in the core, bleed, and heat exchanger are constants but have a 
different (iterated) value in each axial section. 

(4) The many core, bleed, and heat-exchanger gas flow passages are each repre- 
sented by a single equivalent passage. 

(5) Gamma and neutron heating in the hydrogen are negligible. 

(6) Gas axial pressure variation in the core and heat exchanger is linear. 

(7) The thrust and bleed nozzles and the bleed-pump-turbine bypass valve are always 
choked. The topping-turbine bypass valve is never choked. 

(8) The total temperature is constant across each valve. 

(9) The thrust-nozzle chamber temperature equals the core-gas exit temperature. 

(10) The nozzle coolant -tube wall temperature is proportional to the core-gas exit 
temperature. 

(11) The core -gas inlet temperature is proportional to the heat -exchanger exit tem- 
perature. 

(12) Water-pump and bleed-pump turbine exit temperatures are proportional to the 
bleed-gas exit temperature. 
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(13) Heat -transfer from the radiation shield is proportional to the bleed-gas exit 
temperature. 

(14) The bleed gas makes only one pass through the core. 

(15) Inlet conditions to the bleed system are the same as to the main fuel elements. 

(16) At each set of bleed-valves positions, turbine temperature ratios are assigned 
the constant values obtained from the corresponding reference point. 

(17) Hydrogen flow through the nozzles is isentropic. 

Fuel Assembly: 

(1) One equivalent fuel cylinder represents the assembly of concentric fuel rings and 
annular coolant passages in both the main core and bleed system. 

(2) The neutron kinetics are space -independent and are simulated by a point model. 

(3) Neutron kinetics in the bleed system is not simulated. The power generated in 
the bleed fuel elements is assumed proportional to the power produced by the main fuel 
elements. 

(4) Reactivity feedbacks are represented as functions of weighted average core -gas 
and core-water temperatures. 

(5) The .contribution to reactivity of the hydrogen density is negligible. 

Water Flow System: 

(1) The temperature rise across the water pump is negligible. 

(2) Heat transfer in the water plenums is negligible. 

(3) The many water-flow passages in the reflector are simulated by a single equiva- 
lent passage. 

(4) The fractions of total water flow rate passing through the side reflector and both 
inside and outside the flow -divider tubes are constants and equal to the 100-percent de- 
sign values. 

(5) Heat transfer from the bleed-assembly aluminum pressure tubes to the water is 
negligible. 

Heat Exchanger: 

(1) A single equivalent tube with surrounding water flow represents the shell-and- 
tube heat exchanger. 

(2) Neutron and gamma heating in the heat-exchanger tubes and gas are neglected. 
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APPENDIX C 


BASIC EQUATIONS 

The three fundamental equations of fluid flow (continuity, momentum, and energy) 
are reviewed herein and presented in forms suitable for analog computation. 


Continuity Equation 

By considering a mass balance in a stationary volume element through which fluid is 
flowing, the author of reference 7 derived the familiar continuity equation 

— = -V • (pu) (Cl) 

at 

In the present one -dimensional analysis, the form used is 

= - _L (pu) (C2) 

at ax 


Momentum Equation 

Newton's law of motion is applied in reference 7 to a fluid volume: 

a(pu) _ _ g 9P a(pu 2 ) 
at ax ax 


(C3) 


or 


a(pu) = _ g ap + (pu£ dp _ 1 a(pu) 2 
at 3x 2 3x p 3x 

r 9 

o 

Preliminary analog transient runs revealed that the term 3(pu) /3x is negligible. De- 
leting this term and including additional terms to represent friction and sudden expansion 
and contraction yield 
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9(pu) = _ g 9P + (pu) 2 9 P _ 2f(pu) 2 _ K (pu) 2 

3t 9x 2 3x pD p 

r 


(C4) 


Energy Equation 

The energy equation is derived by applying an energy balance to a control volume, 
such as that shown in figure 11(a): 

Change of internal energy in volume = Net heat conducted in 

+ Net heat convected in + Heat generated within 
+ Net work done by pressure forces + Viscous dissipation 



(a) Application of energy balance. 



x 

(b) Application of pressure forces. 
Figure 11. - Control volume. 


In this investigation, conduction within a 
flowing gas is neglected, as is viscous dissipa- 
tion. The small magnitudes of these terms does 
not warrant the amount of analog equipment 
needed to simulate them. 

The total internal energy per unit mass e 
is 


e = E + — (C5) 

2gJ 

where E is the thermodyanmic internal energy 

2 

per unit mass, and u /2gJ is the kinetic energy 
per unit mass. The net inflow of total internal 
energy occurring through the x-faces is 

- JL (pue)dx dy dz dt (C6) 

dx 

in time increment dt. 

The net increase in internal energy within 
the control volume in time dt is 


(C7) 
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— (ep)dx dy dz dt 

at 




Consideration of the pressure forces on the x -faces, as in figure 11(b), leads to an 
expression for the net work done on the volume: 

- A _L (Pu)dx dy dz dt (C8) 

J 3x 


If S is used to denote the sum of the heat-generation rate and heat -convection terms, 
the energy equation can be expressed as 


0(p ue) 3(Pu) t S _ 3 (ep) 

3x J3x V at 


(C9) 


With the use of equations (C2), (C3), and (C5), equation (C9) can be written in the form 


P 



+ 



P Su 
J dx 


In the notation of the total derivative 


D_ a _3_ dx _3_ 
Dt 3x dt at 


equation (CIO) becomes 


DE 
P — 
Dt 


S _ P 0u 
V J 0x 


(CIO) 


(Cll) 


(C12) 


Equation (C12) applies to flow in a constant area passage. If the flow area varies, the 
counterpart of equation (C12) is 


oA 51 = £ _ P 3(uA) 
Dt V J 0x 


(C13) 


as is shown in reference 7. 
For an ideal gas 


E = J ' I Cy dT + const 

= CyT + const (C14) 
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where Cy is assumed constant. Therefore, from equations (C12) and (C14), 

p c v 2T = i-£ 

v Dt V J dx 

With the use of the equation of state for a perfect gas, 

P = pRT 

and the continuity equation (eq. (C2)), equation (C15) can be written as 

p Cp 51 - i =i 

F Dt J Dt V 


(C15) 


(C16) 


(C 17) 


Checkout transient analog runs revealed that omission of the term DP/Dt had a negligible 
effect. Therefore, when equation (C17) is expanded by equation (Cll), it becomes 


3T _ _ u 9T + S 
at dx pCpV 


(C18) 
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APPENDIX D 


APPLICATION OF BASIC EQUATIONS 

In this appendix the fundamental equations are applied to the components used in the 
model. All equations programed in the simulation are presented herein in the form used. 
Stations and nomenclature are indicated in figure 1 (p. 4). Hydrogen properties used in 
the analysis were taken from reference 8. 


HEAT EXCHANGER 
Water Energy 

Averaging . - Of several averaging techniques studied, the following one was con- 
sidered most satisfactory. The steady- state version of equation (C18) was used in aver- 
aging and was applied between the center of a section and its exit. Hence, this method, in 
effect, replaces the three-section model by a six-section model. Use of the steady-state 
version of the energy equation allows solution of the section-exit temperature in terms of 
the section -center temperature. The latter temperature is evaluated by the computer 
from the dynamic differential energy equation. 

Heat -transfer coefficient . - The heat -transfer coefficient used on the water side is 
the familiar Colburn equation: 



This expression is separated into three factors: 


^.(GHM 0 - 8 ^ 


(Dl) 


(D2) 


in which 


0.021 / p f\ 
P 0,2 A 0,8 Ifb/ 


(D3) 


is a constant for each section, and 
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(D4) 


»* = 



varies with the waterside film temperature in each section. A plot of the function 
with film temperature is well approximated by a straight line in the temperature range 
encountered in the heat exchanger. The linear relation derived is expressed as 

= 3. 739X10' 4 T x f - 0. 1247 (D5) 

• • 0 8 

A plot of M against M suggested the use of a linear relation between these two 
variables. In the range applicable to the heat -exchanger water (350 to 1200 lb mass/sec), 
the relation used was 


M®' 8 = 0.212 M y + 38 

X X 


(D6) 


Within this range the maximum error is less than 5 percent. 

Section A. - Application of equation (C18) to section A yields 


( C pf V )x,A%^ = -[“ C P,a( T 2' - T l')] x + [ A S, A G A H A“°' 8 ( T S - T )a] x < D7 > 


which includes the use of equation (D2). With the assumption 


T X, A = 7 ( T l' + T 2') 


at steady state, equation (D7) becomes 



(D8) 
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This expression for T is used in the dynamic version of equation (C18): 

Xj c* 


dT 


M x ,a 


t a + 


^ GA S H \ /m°‘ 8> 


2C, 


M 


( T s - t )a - T r 


— 1 X 


[(ga s h ) a m 0 ' 8 (t s - t ) a ] 3 


When like terms are combined, this equation becomes 


dT x,A_ 

M / T 

- T*^ 

+ 

/ GA s H \ 

dt 

LipVIa 1 

Aj 

X 

_\ 2 c p pv /a J 


-'x 


(D9) 


Section B. - The energy equation applied to section B is 


(<V V )x, B ^ - - [“ C P, b( T 3' - T 2')] x + [(V^ 0 ' 8 ( T S - T )b] x < D10 > 


Solution of this equation at steady state for T x g» gives 


T x, 3' " 


T B + 


f'GA c H\ /,>0. 8\ . x 

4 r) (Vfe • T > 


B 


-Jx 


(Dll) 


This relation and that derived in equation (D8) for T x ^ are substituted into equation 
(DIO): 


dT. 


x, B _ 
dt 


M (T A - TJ + ' GAgH 


(PV) B ’ A B \2CpPV 


M 


0.8 


( T s - T )j 


'GAgH 
l2CppV, 


M 


0.8 


( T s - T ) 


B 


B 


— *x 


(D12) 
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Section C. - In section C, the energy equation is 


dT x £ 

Cp pV c — = 

P X , L dt 


= [-“ C P, c( T E ' T 3')] x + [f A S GH )c “°' 8 ( T S - T )c] x < D13 > 


The heat-exchanger water-exit temperature is evaluated from the steady-state version: 


l x,E 


T C + 


/GAgH' 

2CI 


!M 


0.8N 


M 


( T s - T )< 


CD 14) 


This equation together with (Dll) is used in equation (D13) so that 


dT. 


x, C _ 
dt 


M 


(PV), 


(t b -t c ) + 


'GAgH 
\2C^pVj 


M 


0.8 


( T S ‘ T ) 


B 


'GAgH 

2CpPV i 


M 


0.8 


( T s - T )< 


(D15) 


Gas Energy 

Averaging. - The averaging technique used is the same as that used for the heat- 
exchanger water. 

Heat -transfer coefficient. - The convective heat-transfer coefficient used to express 
the transfer of heat from wall to hydrogen was obtained from a correlation presented in 
reference 9: 


h X = ' 


0.021 
D °. 2 a 0 . 8 




0. 29+0. 0019(x/D) 



This equation is separated into the following factors: 

h X = (HbM°- 8 )x 


(D16) 


(D17) 
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where 


and 



(D18) 




29+0. 0019(x/D) x 


(D19) 


The length x in the last equation is the distance from the heat-exchanger inlet to the 
center of a section. The factors H x and b x were each assigned a constant value that 
differed for each section. Several iterations were made by using gas and wall tempera- 
tures determined from analog runs before final values of the constants were specified. 

. o Q 

The linear relation used for the factor M x was 

m£ 8 = 0. 34 M x + 6 (D20) 

This approximation is close to an exact value at hydrogen mass-flow rates between 50 
and 105 percent of design. At 36. 6 percent (the bleed-pump stall limit), it is 4. 7 percent 
too high. 

The form of the energy equation used for the heat -exchanger gas is 


dTx 

dt 


RT X 

V X P X 


hA„ 

- M AT + — ^ (T q - T) 
C & 

^■n 


X 


(D21) 


It includes the use of the perfect-gas equation of state (eq. (C16)). Application of the 
averaging method described in the heat-exchanger water-energy analysis yields the 
following equations for the three sections: 
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Section B: 



Section C: 



In the course of the derivation of equation (D24), the steady-state form of the energy 
equation yielded 


L X, E 



-■x 


(D25) 


Gas-pressure distribution. - Because only one axial section was used for the heat- 
exchanger gas momentum equation, a linear pressure variation was assumed to evaluate 
the pressure at the center of a section. The following relations were used: 


AP ■ ( p x ' p cn) e 

(D26) 

P X, A = P CN, E + ^ AP X = \ P CN, E + g P X, E 

(D27) 

P X, B = P CN, E AP X = \ ( P CN + P x) E 

(D28) 

P X, C = P CN, E + ^ AP X = (g P CN + ^ 

(D29) 
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Gas Momentum 


For the one axial section employed to simulate the momentum equation applied to the 
heat -exchanger gas, the form of equation (C4) used was 


dM } 

dt 


X,S 

Ax 


/He 


AP X + I-1 

x \PV, 


' Ap ' 


f2f Ax) 


'X 


He 


(D30) 


Calculation of the following parameters is necessary for insertion into equation (D30): 


Ap x - (p x - P C n) e 

(D31) 

p X = \ ( p CN + p x) E 

(D32) 

' X > E - 1 (M x - M tt - Mg) 

dt v x _ T 

(D33) 

AP x = ( p x " p cn) e 

(D34) 

P X, E = p x, E RT X, E 

(D35) 

f _ 0.84Z 

(D36) 


(In Re x f - 1. II) 2 ' 22 


In equation (D36), the factor Z represents the friction entrance and exit effects due to 
contraction as the hydrogen gas enters the heat -exchanger tubes and expansion as it 
emerges. The remainder of the expression in equation (D36) was derived from a plot of 
the Karman-Nikuradse approximation for friction in fluid flow. 


Tube 

Averaging. - Within the tube material, the radial temperature variation was assumed 
to be linear. With representing the temperature at the center of the section on the 

water side and T R the center temperature on the cold hydrogen side, the surface tem- 
peratures are expressed as follows: 
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Waterside surface: 


T x , s ® T w " ® Tpj (D37) 

Gas side surface: 

t x,s =1 - 5t h- 0 5 t w (e>38) 

These equations were used in the heat -exchanger water analysis (eqs. (D9), (D12), (D14), 
and (D15)) and gas analysis (eqs. (D22) to (D25)). 

Energy equation . - A heat balance applied to the waterside section yields 

( pvc p)„ e ? ■ V*, s( T - T s) x - (y^w - t h> 

and to the gas side is 

( pvc p) He = _h x A x, s( T s ' T ) x + (T w " t h) (D40 > 


Application of these two equations to axial section A gives 



The equations for axial sections B and C are identical to the preceding two equations ex- 
cept that subscript A is replaced by B and C, respectively. 


60 



CORE 


Gas Energy 


Heat -transfer coefficient . - The form of equation (D16) was used again to express 
the convective heat -transfer coefficient between the hydrogen gas and fuel surface. In 
this instance, it was separated into two factors: 


h 


c 



The factor 



(D43) 


(D44) 


is a constant that has a different value for each of the three sections. Several iterations 

were made to establish the final values used. Equation (D20) was used again to evaluate 
• 0 ft 

M^' and has the same range of accuracy. 

Section A . - The basic form of the energy equation used is similar to that applied to 
the heat- exchanger gas. 



When dT 
sion for 


c ^/dt is set equal to zero in equation 
T c 2 ’ in to this equation yields 


(D45), insertion of the resultant expres- 



The value of 



is determined by 



(D47) 
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The approximation in equation (D47) is good because the small temperature drop across 
the topping turbine subtracts from the small hydrogen gas temperature rise in the shield. * 
The factor (T t /T y tp) > which is slightly greater than 1, was assigned a different 

C, 1 A, Hi re f 

reference value for each set of bleed-valves settings. 

Section B. - After the steady-state energy equation is solved for T c g», the dynamic 
version becomes 



Section C- - The corresponding equation for section C is 



The core-gas exit temperature T c g, determined in deriving equation (D49), is 


c, E T c, C + 


/A F, S H c \ l M® 1 8 ' 


2C 


P,c, 


M 


( T F - T c), 


(D50) 


Gas-pressure distribution . - The core-gas-pressure distribution was assumed to be 
linear because only one axial section was used in the mass-flow rate calculation: 


AP c " P c,E _ P c,I 


= P„ T + - AP, 


c, A " c, I 


(D51) 

(D52) 
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(D53) 


P 


c, B 



AP c 


P 


c, C 



+ - AP 
4 C 


(D54) 


Gas Momentum 


Application of the momentum equation (C4) to the one-axial section results in 



(D55) 


The last term represents the entrance and exit losses of the individual axial stages of the 
fuel elements. The following parameters were evaluated in solving equation (D55): 

Ap c = ( P E ' p l) ( D56 ) 

c 


^c, E 
dt 


1 


V c-N 


(M c - M n ) 


(D57) 


Preliminary digital calculations of a multisection model of the core gas revealed that the 
axial variation of 1/p is nearly linear. The relation used is 


p c 



(D58) 


Since two points determine a straight line, l/p c ^ could have been used in this expres- 
sion instead of l/p c B or l/p c c with an appropriate adjustment in the constants. 

The core -gas density factors were evaluated by 



(D59) 


and 
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1 


(D60) 



f 0.84 Z 

C (In Re c - 1. II) 2 - 22 

Core-gas exit pressure was calculated from 

P c, E = Pc, E RT c, E 


(D61) • 


(D62) 


FUEL 


A heat balance in the cylindrical model of the fuel elements included heat generation. 
The three axial sections used corresponded to those in the core-gas analysis. The heat- 
generation rates in each section were based on the axial fuel distribution rates: 

Section A: 


Section B: 



dt 




Section C: 




^F, B “ °* 325 ^TOT 



Qp q — 0. 365 Q'pQ'p 


(D63) 

(D64) 

(D65) 

(D66) 

(D67) 

(D68) 
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MODERATOR WATER 


The moderator water flows through the core in two annular passages separated by an 
aluminum flow divider. To distinguish between the two moderator waters, the volumes 
inside and outside the flow divider are referred to herein as "inside” and "outside" water. 


Heat-Transfer Coefficient 


The convective heat -transfer coefficient used was the same as that used for the heat- 
exchanger water: 




* 

i 0.021 

feV'V 6 

/n \ o. jT 
m0. 8/ f\ l 

d 0.2 a 0.8 

w J 

f \ p b/ 

' f jy 


\, = (™°' 8 )w 


(D69) 


(D70) 


The factor 


H w = 


0.021 
ID 0, 2 A 0 * 8 



0.4 


, 0.6 



(D71) 


was assigned a constant value, determined after several iterations, for each passage. 

For the range of flow rate of both the inside and outside water, the following linear 
• 0 8 

expressions for M ‘ were used: 


<* = 0. 235 25 

<o=»-235 M^o+aSJ 


(D72) 


This relation is applicable in the flow rate range between 200 and 800 pounds mass per 
second, within which the maximum error is 3. 6 percent. Throughout the analysis, the 
inside and outside water mass-flow rates were assumed to be constant fractions of the 
total, the constants being evaluated at the 100 -percent design point: 
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(D73) 


M w, i = °' 423 M w 


M w>o = 0 - 4 87 M w 


(D74) 


The two heat-transfer coefficients on each of the inside water boundaries were assumed 
equal. 


h D, i “ h PT, o " h i 


(D75) 


Energy Equation for Inside Water 

Heat enters the inside water by convection from the pressure-tube and flow -divider 
surfaces and by internal generation: 


dT. 


^PT.sVI mO.S, T wi ) 

* \P V / w,t KpL , ’ ’ 

W, 1 




' W Jx 


(pV Cp ) 

AT w, i “ 2 < T w, i " T x, 
Qwj i = °* 00598 Q'pQ'p 


w,i u D w,i' + l nVf 


pVC, 


w,i 


(D76) 

(D77) 

(D78) 


Energy Equation for Outside Water 

Only one surface, the flow-divider outer surface, supplies convected heat to the out- 
side water. Therefore, 


dT w,o. I 

(m AT^ 

\ + 
W, 0 

A D,S«w 

dt ' 

Vpv ) 

1 

KpI 




-JO 


w, o 


AT w, o 2 < T w, o ” T x, 

Qw,o = 0- 01600 Qtot 


(D79) 

(D80) 

(D81) 
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PRESSURE TUBE 


Heat enters the pressure tube by conduction through the stagnant hydrogen layer into 
. its inner surface and by heat generation within its material. Heat leaves the pressure 
tube outer surface by convection into the inside water. A heat balance applied to the 
pressure tube gives 


dT pT _ A pT> gH^ i 

lvl w, r A w, i 


dt 


(pVCp) 


PT 


) + 


Q sh + Q 


PT 


(pVCp) 


PT 


PT 


(D82) 


The heat loss from the core^gas conducted through the stagnant hydrogen layer Qg H is 
assumed to be a linear function of the total heat rate produced. It also includes the effect 
of thermal radiation. This heat loss was evaluated at 20 and 100 percent to determine the 
function: 


Q sh = 0. 00255 Q tqt + 2879 (D83) 

The heat generated within the aluminum material of the pressure tube is proportional to 
the total heat rate: 


Qprp — 0.00218 Q'pQ'j’ 


(D84) 


FLOW DIVIDER 


Heat is generated within the material and leaves through its inner and outer surfaces 
by convection into the inside and outside water, respectively. The heat balance is 


dT~ ( A D,S H w). n q ( A D,S H w) n o / n 

D _ ’ 1 M°- 8 AT . - T ) + 0 M°- 8 (T - T.0 + I Q 

‘^■’117 \ ' ^ TIT 1 ' T>> + /„,rr \ W,0 V W, O" D' 


dt 


(PVC P ) D (PVC P ) D 


pVC, 


(D85) 


'D 


Qp — 0. 00132 QrpQrp 


(D86) 
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SIDE REFLECTOR 


Heat-Transfer Coefficient 

The convective heat -transfer coefficient used is defined by equations (D69), (D70), 

• 0 8 

and (D71) with the subscript R preceding w. The linear relation used for M R ^ is 

m"; 8 w =°. 3 4 M RjW + 6 (D87) 

where 

^R,w = 0 - 090 k w (D88) 

The linear expression in equation (D87) is applicable in the mass-flow rate range of 35 to 
120 pounds mass per second, within which the maximum deviation is 3. 5 percent. 


Water Energy 


Heat enters the reflector water by convection from the reflector wall and by genera- 
tion within the water. The energy equation is 


dT 


R, w _ 
dt 


[ M AT 

V pv 


w 


A S»w 

( PVC P) W 


M w 8 < T - V + 


Q 


[pv C, 


w 


(D89) 


-*R 


^R, w 2 ( T R,w- T x, E ) 

(D90) 

^R, w = 4 * 57x10 4 Q-tot 

(D91) 


Reflector Material 


Heat is generated within the reflector material and is convected from its wall into 
the reflector water. The heat balance is 

^R 

dt 

Q r = 0. 00340 Q-tot (D93) 


A . sHw M°- 8 (T - T) + 

PVC D w 


Q_ 

lPVC n 


(D92) 
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REFLECTOR- AND CORE-WATER EXIT MIXING 


The three water streams of the reactor mix in an outlet plenum. For a heat balance 
.in the plenum, 


dT 1 /, . . . Qy 

^w, i^w, i, E + ^w, o^w, o, E + ^R, w^Tt, w, E ” Vmmw + 


dt (pV) 


mmw 


'P, mmw/ 
(D94) 


The term is the heat generated in the heat -exchanger water; it is lumped with the 
water-plenum mixing for simplicity. This parameter is evaluated as 


-4 


Q* = 4. 24X10 Q tqt 


(D95) 


WATER-FLOW LOOP 

Application of the momentum equation to the water-flow loop yields 

dM a 2 

_Z = IA_ (AP wp - AP f ) (D96) 

dt V ’ 

w 

The friction pressure drop in the water loop is assumed to be 

AP w>f r= °. 02037 (D97) 

BLEED SYSTEM 
Gas Energy 

The heat -transfer coefficient of reference 9 was also used in the bleed system analy- 
sis. It is expressed as 


h BG " 



0.6 


J BG 



0. 29+0.0019(x/D) bG 



(D98) 
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and 


h BG = (“M 0 ' 8 )bg < D99 > 

In the bleed system, H BG , which is a function of bleed fuel geometry, gas properties, 

and bleed-gas temperature, was assumed constant. Several iterations preceded the final 

• 0 8 

determination of the constant. A plot of M BG against Mg G in the range between 0. 8 
and 3. 3 pounds mass per second, as used in the bleed system, is essentially linear. The 
form used is 


= 0. 0235 ]fcl BG + 0. 3133 


(D100) 


The maximum error within this range is 5 percent. 

Equation (C18), when applied to the one-section model of the bleed system, yields 


dT 


BG 


dt 


= E (pv) bg 


■<M AT) bg 


B, S H BG\v,0. 8 


P,BG 


m bg< t bf - T ^-> 


BG' 


(D101) 


In this equation 


AT BG = T B, E " T c, I 


(D102) 


and 


’bG “ ~ ( p c, I + P B, E^ 


(D103) 


Gas Momentum 

The momentum equation (C4) applied to the bleed system is 



The following relations are used in this expression: 


(D104) 
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AP BG " P B, E ' P c, I 


Ap 


BG _P 1 " p c,I 


P B,E = P 1 + 0-1097 MbG 


P BG " \ ^1 + p c, 1^ 


BG 


0.84 Z 


ln( Re B G,f ' 1 * 11 ) 


2.22 


Fuel 


A heat balance applied to the solid fuel materials yields 


(D105) 

(D106) 

(D107) 

(D108) 

(D109) 


dT BF 

dt 



^BF “ 0*0101 Q-rOT 


(duo) 


(Dill) 


The approximation of equation (Dill) was used to eliminate the necessity of simulating 
the kinetics equations in the bleed system. 


NOZZLES 

Nozzle Coolant 


Energy equation . - The following equation was applied to calculate the gas tempera- 
ture at the center of the nozzle-coolant section. 



v. 




(D112) 
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In this equation, the convective heat-transfer coefficient is expressed as 


h CN = 0.0182 M°-J 

1 

(D113) 

where 

* 

Mcn = °- 34 m C n + 6 

(D114) 

as used in the heat-exchanger and core-gas simulations, and 


df> CN,E = / 1 Vfa -fa ) 

* \ V CN-x/ 

(D115) 

P CN = \ ( p I + P e) cn 

(D116) 

T CNW = °- 2108 T c , E 

(D117) 

T CN, E = 2T CN " T CN, I 

(D118) 

The nozzle -coolant inlet temperature T^ N j was assigned a different constant value for 


each set of bleed- valves settings. 

Momentum equation. - Application of equation (C4) yields 


dM 


CN 


dt 


= S 


AP’ 


3. 486 fa 


VAX 


' CN 


CN 


(D119) 


in which 


-AP cn - P„ + AP DT1 + AP 


BP 


TP 


CN, E 


(D120) 


and 


P CN, E ~ P CN, E RT CN, E 


(D121) 
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Nozzle-Gas Flow 


Both the thrust and bleed nozzles are assumed to be choked throughout the analysis. 
For isentropic flow the mass-flow rate is 


M 


-a P Jyg/ 2 \r*Vr-l 

! R \r + V 


(D122) 


With the use of the perfect gas relation, this equation becomes for the thrust nozzle 


'*„* 107 - 7 Vn,lfl (D123) 

and for the bleed nozzle 

Sd BN = 107. 7 A BN p 4>E yi- (D124) 


PUMPS 


A conventional pump-performance map relating pressure rise and flow rate can be 

2 

compressed into a single curve when plotted as AP/N against M/N. Likewise, one 
2 

curve relates L/N with M/N. For the three pumps, the six variations on such plots 
were set on variable diode function generators, as indicated in the following equations: 


Bleed pump: 



(D125) 


(D126) 
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Topping -driven pump- 


Water pump: 



(D127) 


(D128) 


(D129) 


(D130) 


TURBINES 


Water- Pump Turbine 

Torque . - The relation between L^^/M^rp^T^ and N^rp/^T^ is linear for a 

constant pressure ratio. In addition, a plot of L^^/M^rp^T^ against the pressure- 
ratio function 



is linear for the zero-speed condition. When these two linear relations are combined, 
the following equations result: 
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or 


Lwt " 


= M, 


WT 


16.46 


f WT 


- 0. 0045 N, 


WT 


(D131) 


In equation (D131), the function is provided by a variable diode function generator 



Mass -flow rate. - The function 



(D132) 


(D133) 


was set on a variable diode function generator. Water-pump turbine mass-flow rate was 
determined from 



(D134) 


Speed. - The turbine and pump torques are used in calculating rotational speed: 


dN- 


WT 


30 


dt 


1 WTpy 


(Lwt ■ 


(D135) 


Hydrogen Bleed-Pump Turbine 

Torque . - The torque -parameter plot of the hydrogen bleed-pump turbine also has 
linear characteristics, in the manner discussed in connection with the water-pump tur- 
bine. The linear equation obtained is 




(D136) 
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Again, a variable diode function generator was set with the same pressure ratio function 





(D137). 


as that used for the water-pump turbine. Accordingly, the form of the torque equation 
used is 


L BT ~ ^Bt( 22- 1 V T 3,E f BT " °' 0082 N Bt) 
The value of Tg g used was approximated by 



(D138) 


(D139) 


in which (T, E /T, ) was assigned the design value when the bleed valves were set at 

the 100 percent positions. The values at the other two reference points were used for 
the 60 and 36. 6 percent settings. The temperature drop across the water-pump turbine 
is less than 4 percent; thus, the assumption in equation (D139) is very good. The pres- 
sure Pg g, used in equation (D136), is determined from equation (D170). The pressure 
P4 g is calculated by 

P 4, E = p 4, E RT 4, E (D140) 

in which 



(D141) 


and 


»4,E 

dt 


(^BT + 


4-N 


BN 7 


(D142) 


The assumption of constant temperature ratio across the bleed-pump turbine (eq. (D141)) 
is not as valid as that used in the water-pump turbine because the former turbine has a 
greater work output. At the 100 percent design point, the temperature at the bleed-pump 
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turbine inlet drops by 14 percent. Smaller percentage drops occur at lower power levels. 

.The 14 percent is not an error, of course. Furthermore, a different reference value is 

used for (t 4 p/T.) at each of the three bleed-valves settings. 
v ’ v ref 

» Mass-flow rate . - Because the bleed-pump turbine is choked at design conditions, 
the following relation was used: 


m BT 


0.00386 


P 3,E 

V T 3,E 


With the use of the equation of state (eq. (C16)), this equation becomes 


(D143) 


“bt= 296 P3,eVVe < d144 > 

This approximation reduces the amoimt of analog equipment required. A 2-percent error 
occurs if the turbine pressure ratio drops from 5. 5 to 4. 0. An error of this magnitude 
is within the accuracy of analog work. 

Speed . - The dynamic equation for turbine (and pump) speed is similar to that used 
for the water-pump turbine: 


pp = -p- < l bt - L BP> < D145 ) 

at ^BTP 

Topping Turbine 

Torque . - The following linear relation also proved applicable to the topping turbine 
torque representation: 


L TT 

M T tV T X, E 



A Nrprp 

9.36X10 ^ — L A - 

V t x,e 


(D146) 


As for the other turbines, the pressure -ratio function was set on a variable diode func- 
tion generator 



(D147) 
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Thus, the topping-turbine torque was simulated by 


* 

L tt = M TT ( 8 . 8 yT X E fp T - 9.36X10" 4 N TT ) (D148) 

In this equation it is assumed that 


T TT, I " T X, E 


(D149) 

Mass-flow rate. - A variable diode function generator was set with the function 


,» _ M T tV T X, e 


TT 


TT, I 


The mass-flow rate was calculated from 


m xt = -^== f tt 


rx, 


E 


Speed. - The same rotational speed equation is used 


dN TT 

dt 


30 


TTPy 


(Ltt - Lipp) 


(D150) 


(D151) 


Exit pressure. - Topping -turbine exit pressure was calculated by adding a friction 
pressure change between core inlet and turbine exit. This pressure drop was assumed to 
be proportional to the square of the total hydrogen mass-flow rate. 

P TT, E = P c, I + 1 * 04 (m tt + IVI 5) 2 (D152) 

A different value of the constant in this equation was used at each of the three bleed- 
valves settings. 

The core-gas inlet pressure was calculated by using the following relations: 


dp 


M = 

dt V 


(Mpp + Mg - M c - Mjjq) 


T-c 


(D153) 

(D154) 
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p c,I = ‘>c,I RT c,I 
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Figure 12. - Typical experimental butterfly- 
valve characteristics. 


The range of pressures and temperatures and 
the absence of a positive shutoff requirement permit- 
ted the use of butterfly -type valves. The need of se- 
lecting a valve size was eliminated by employing a 
relation between valve flow -capacity coefficient ratio 
c v / c v max ancl stem position. A plot of the charac- 
teristics of a typical experimental butterfly valve is 
presented in figure 12. 

The following equations were used to calculate 
the mass-flow rate passed by a valve: 


For nonchoked flow: 


M 


0.008 



For choked flow: 


(D155) 


M = 0.004 



(D156) 


In the application of these equations to the computer, a relay was used with the simula- 
tion of the two water -pump turbine valves and the bleed-pump-turbine throttle valve to 
provide for choking or nonchoking flow. Valve pressure ratio was used as a basis of the 
relay operation. In the following explicit mass-flow rate equations for these valves, only 
the nonchoked relation is given, it being understood that the relay makes the appropriate 
selection. Because the bleed-turbine bypass valve is choked at the three reference 
points, it is assumed to be choked throughout the simulation. Also, the topping-turbine 
bypass valve is assumed to be nonchoked throughout because it is not choked at the design 
point. 


Water-Pump-Turbine Bypass Valve 

This valve is designated herein as valve 1. The mass -flow rate it passes is 


79 



Mj = 0.008 


P 

yl t 



The following variables are calculated for use in this equation: 



1 


V 


B-l 




-m 2 ) 


T 


1 " 



T 


BG 


(D157) 


(D158) 


(D159) 


This equation represents by a constant temperature ratio the temperature drop of the 
bleed gas as it passes through the shield. The same constant applies to all three refer- 
ence points: 


P 1 = P 1 RT 1 


dt 


(M^ - M^yrp) 


1-WT 


P 


1,E ” p l, E RT 1 


Water- Pump-Turbine Throttle Valve 

The mass-flow rate is expressed as 


M 2 = 0.008 c y 2 





The thermodynamic variables required in this equation are obtained from 


(D160) 

(D161) 


(D162) 


(D163) 
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(D164) 


Po = — P, 


'ref 


dp 


2, E _ 1 


dt V 


2-3 


O^wt + Mg - Mg - M 4 ) 


T 2,E “ T 3,E 


E “ p 2, E RT 2, E 


Bleed-Pump-Turbine Throttle Valve 


The mass -flow rate is 


M 3 = 0.008 c v 3 


P 2 *£3 _/^3 

f P » 


-<E 


dp 


3, E 1_ 

dt V 


3-BT 


(M^ " M-grp) 


P 3, E = p 3,E RT 2,E 


Bleed-Pump-Turbine Bypass Valve 


The choked equation applied to this valve is 


M 4 = 0. 004 c y> 4 




and 


P 4 = 


lP / 2,E 

L 2. E /ref 


(D165) 

(D166) 

(D167) 


(D168) 

(D169) 

(D170) 


(D171) 


(D172) 
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Topping-Turbine Bypass Valve 


* 

For this valve, the nonchoked equation is 


M 5 


0. 008 c y 5 — T ’ -■ 


f*X, 


E 


P TT, I " P X, E 



(D173) 


(D174) 


The value of Krprp was changed to match reference conditions at each of the three sets 
of bleed- valves positions. The value for P,p^, is obtained from equation (D152). 


REACTOR KINETICS 

The neutron kinetics is represented by a one-energy group model with six groups of 
delayed neutrons: 


dn _ 
dt 



where 


6k' = 5kp + 6k^ + 6kp 


(D175) 


(D176) 


The six equations describing the net rate of formation of the delayed neutron precursors 
are of the form 


dC., n/3,, 

— = -4- - X V C V i’ = 1, 2, ... 6 (D177) 

dt l* 11 

The variables in equations (D175) to (D177) are assumed to be independent of space. 

A distinctive characteristic of the reactor is the negative temperature reactivity fed 
back by both the fuel and the water temperatures. The positive contribution to reactivity 
of the hydrogen density in the core is negligible by comparison. 

In the fuel-tungsten mixture a rise in temperature causes the resonance region to be 
broadened; this phenomenon is known as the neutron Doppler effect. Therefore, more 
neutrons are absorbed at high energy, which leaves fewer to be thermalized for subse- 
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quent fission and, thus, reduces the reactivity. The reactivity fed back by the average 
fuel temperature was expressed as 

5kp = -2. 67X10" 6 T p - 0. 01335 (D178) 

The average fuel temperature was approximated by the temperatures of the three axial 
fuel sections, each weighted approximately according to the reactor-power distribution 
in its section; that is, 

T f = 0. 299 T Fj A + 0. 339 T F B + 0. 362 T p c (D179) 

The water-temperature reactivity could not be approximated well by a linear function 
of temperature. Therefore, the variation was set on a variable diode function generator: 

6k^ = f”(T w ) (D180) 

The average water temperature T w was determined from the moderator inside and out- 
side water temperatures. The weighting contribution of each was assumed to be propor- 
tional to its volume: 

T w = 0. 2734 T Wj . + 0. 7226 T w> Q (D181) 
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